A Theory for the High-T c Cuprates: Anomalous Normal-State and 
Spectroscopic Properties, Phase Diagram, and Pairing 
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A theory of highly correlated layered superconducting materials is applied for the cuprates. Dif- 
ferently from an independent-electron approximation, their low-energy excitations are approached in 
terms of auxiliary particles representing combinations of atomic-like electron configurations, where 
the introduction of a Lagrange Bose field enables treating them as bosons or fermions. The energy 
spectrum of this field accounts for the tendency of hole-doped cuprates to form stripe-like inhomo- 
geneities. Consequently, it induces a different analytical behavior for auxiliary particles correspond- 
ing to "antinodal" and "nodal" electrons, enabling the existence of different pairing temperatures 
at T* and T c . This theory correctly describes the observed phase diagram of the cuprates, including 
the non-Fermi-liquid to FL crossover in the normal state, the existence of Fermi arcs below T* and 
of a "marginal- FL" critical behavior above it. The qualitative anomalous behavior of numerous 
physical quantities is accounted for, including kink- and waterfall-like spectral features, the drop in 
the scattering rates below T* and more radically below T c , and an effective increase in the density 
of carriers with T and uj, reflected in transport, optical and other properties. Also is explained the 
correspondence between T c , the resonance-mode energy, and the "nodal gap". 

Keywords: superconductivity, cuprates, auxiliary particles, anomalies, pairing 



I. INTRODUCTION 

High-T c superconductivity (SC) has been in the fore- 
front of condensed-matter physics research since its dis- 
covery in the cuprates over 23 years ago. This system is 
characterized by anomalous behavior of many of its phys- 
ical properties [MV\\ which led to the suggestion of non- 
Fermi-liquid (non-FL) models fl2T - [l5j for its electronic 
structure. The recent discovery of high-T c SC in iron- 
based compounds (referred to below as FeSCs), provides 
a new test case for high-T c theories, especially in view [l6[ 
of the striking similarity of their anomalous properties to 
those of the cuprates. 

Recently, a unified theory for the cuprates and the 
FeSCs has been derived [16| on the basis of common 
features in their electronic structures, including quasi- 
two-dimensionality, and the large- U nature of the elec- 
tron orbitals close to the Fermi level (E w ). Electrons 
of such orbitals cannot be treated through a mean-field 
independent-particle approach; therefore, their behavior 
is studied in terms of auxiliary particles, representing 
combinations of atomic-like electron configurations. 

Within this auxiliary-particles method [17| , configura- 
tions in sites i can be, approximately, treated as bosons 
or fermions, if each site is mathematically constrained to 
be occupied by one (and only one) configuration. There 
is a freedom in choosing configurations of an odd num- 
ber of electrons as fermions, and of an even number of 
electrons as bosons, or vice versa [161 ]. The choice here 
is that configurations corresponding to undoped cuprates 
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or FeSCs (thus the parent compounds) are bosons. 

These auxiliary particles have been treated beyond 
mean-field theory [14]]. A grand-canonical Hamiltonian 
rt has been written down 16], where (in addition to 
the chemical potential fj,) a field of Lagrange multipli- 
ers Xi was introduced (A = (Aj)) to maintain the above 
auxiliary-particles' constraint. This Lagrange field rep- 
resents an effective fluctuating potential which enables 
the treatment of the above configurations as bosons or 
fermions; its effect on them is analogous to the effect of 
vibrating atoms on electrons, and similarly to lattice dy- 
namics, this is a Bose field [Iff (see discussion about it 
below). 

These auxiliary particles are assigned the following 
ad hoc names [16| : (i) combinations of boson atomic- 
like configurations are referred to as "svivons"; their 
Bose condensation results in static or dynamical inhomo- 
geneities which could be manifested in the existence of 
a commensurate or an incommensurate resonance mode 
Hi : {it) combinations of fermion atomic-like configura- 
tions are referred to as "quasi-electrons" (QEs); they 
carry charge — e and are introduced through electron or 
hole doping (in the second case the configurations corre- 
spond to holes of QEs); (Hi) the Lagrange-field bosons 
are referred to as "lagrons", and their coupling to the 
svivons and QEs provides a dynamical enforcement of 
the auxiliary-particles' constraint. 

Stripe-like inhomogeneities have been observed [l8r - 
[22|] within the phase diagram of the cuprates. Even 
though such an inhomogeneous behavior is often dynam- 
ical within the SC regime, it is not treated here as a 
perturbation to a homogeneous state, but as a correction 
to a static inhomogeneous structure. The application of 
the constraint, approached here through the Lagrange 
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field, must take into account the existence of these inho- 
mogeneities. This results in a multi-component scenario. 

In small-t/ systems, where low-energy electrons are ap- 
propriate quasiparticlcs (QPs), the FLEX approximation 
[23l | could be applied [2J| to derive an effective pairing 
interaction between them, due to their own spin and 
charge fluctuations. The low-energy large- U electrons, 
discussed here, cannot be considered as appropriate QPs; 
but their description in terms of independent auxiliary 
Fermi and Bose helds provides an intrinsic pairing mech- 
anism through QE-lagron coupling; it can be treated [iH ] 
within the Migdal-Eliashberg theory [25T - l27j , similarly to 
electron-phonon coupling in the case of phonon-induced 
pairing, but with coupling constants which are so strong 
that lattice instabilities would have been driven in the 
electron-phonon case. 

In this paper auxiliary particles which are specific for 
the cuprates are introduced, and their Green's functions 
are applied to derive those of the electrons. The detailed 
QE and electron spectral functions are calculated, with 
further emphasis on the anomalous low-energy features. 
The QE, svivon and electron scattering rates are derived 
in the non-FL regime, and the resulting physical anoma- 
lies are discussed. 'K is then applied to derive QE and 
electron pairing which is demonstrated to result in the 
distinct pseudogap (PG) and SC phases; the associated 
anomalous physical behavior is discussed. Further calcu- 
lations on part of these aspects (specifically in the paired 
states) will appear in separate papers. 



II. HAMILTONIAN 

Electronic-structure calculations, and a variety of ex- 
perimental data, indicate that the "universal" (thus 
not characteristic of just specific compounds) low-energy 
properties of the cuprates are primarily determined by 
Cu(d) and 0(p) orbitals of the Cu02 planes. Within 
the LDA, these orbitals generate a multiple-band struc- 
ture extending over a range of [28[ about 9 eV. However, 
the electronic states at the vicinity of E F dominantly 
reside I28l - l30t | in one of these bands, of Wannier func- 
tions 31] which are centered at the planar Cu atoms, 
and correspond mainly to er-antibonding states between 
Cu(d x 2 _ y 2) and 0(p x ,p y ) orbitals (though a minor con- 
tribution of other orbitals exists as well). Due to this 
orbital nature of the band, its treatment requires effec- 
tive transfer (hopping) integrals up to the third-nearest 
neighbor (see below). 

The effective intra-site Coulomb parameter U, for elec- 
trons in this band, corresponds to its Wannier functions 
which are, primarily, combinations of Cu(d) and O(p) 
functions. Estimates of U must take into account the fact 
that the d intra-atomic integral Ud is considerably larger 
than the p integral U p . In an early derivation of the effec- 
tive one-band Hamiltonian [32| i n hole-doped cuprates, 
the large-[/,j and small-J7 p limits have been assumed, re- 
sulting in approximately ionic Cu(cZ) and itinerant O(p) 



states. The resulting [32j upper Hubbard band corre- 
sponds to Cu +1 ions, and the lower Hubbard band to 
Cu +2 ions, part of which are forming singlet states with 
the doped 0(p) holes (known as "Zhang-Rice singlets"). 

Consequently, the value of U, derived through this 
analysis, consists of Ud plus corrections, due to the sin- 
glet energy and the difference between the one-electron p 
and d energies which are considered to be substantially 
smaller. A more realistic evaluation of U would result in 
a somewhat smaller value than Ud', however, since the 
width of the effective band, within a two-dimensional 
(2D) Brillouin zone (BZ), is [HrHl about 3 eV, it is 
appropriate to treat it using a large-?/ approach. 

Within this one-band approach, a 2D square lattice of 
N sites i is considered (assuming translational symmetry 
between them); the Cu atoms are located at the lattice 
points R;, generated by the unit vectors ax and ay. An 
electron of spin a in site i is created by d\ (the two 
spin states are denoted as a =t,4-, or a = ±1). The 
Hamiltonian is expressed as: 

where the transfer integrals t(R) are considered up to 
the third-nearest neighbor, and expressed in terms of the 
parameters t, t 1 and t": 

t(±ax) = t(±ay) = -t, (2) 
t(±ax ± ay) = -t', t(±2ax) = t(±2ay) = -t". 

The present auxiliary-particles approach [lj| then be- 
comes the "slave-fermion" method, applied in previous 
works by the author [33 35] . This approach is different 
from homogeneous RVB models jl2l . H . [ll| which are 
also based on auxiliary particles, but generally within 
the "slave-boson" method, where the fermions are spin- 
carrying "spinons" and the bosons are charge-carrying 
"holons". Within the present approach, the roles of 
fermions and bosons are switched, and consequently in- 
homogeneities are inherently introduced through Bose 
condensation. An electron creation operator is then ex- 
pressed as [H, E3] : 

4x = 4a + cre K-<T> ( 3 ) 

where s icr , h\ and ej create a boson svivon state, and 
fermion holon and "doublon" states, respectivel y. The 
auxiliary-particles' constraint is expressed as [l6llr7j: 

E a *W+^+eK = l. (4) 

For hole-doped cuprates, the effect of the upper- 
Hubbard-band doublons is ignored in the lowest order, 
and the creation operator of a fermion QE state is de- 
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fined as: q\ = h i . Eqs. (|3I4|) are then approximated as: 



E4 



(5) 
(6) 



By expressing the electron operators in the rhs of Eq. ([T]) 
in terms of the auxiliary-particle operators, through 
Eq. ([5]), one gets that the U term in H d drops, as is 
expected far below the upper Hubbard band. 

Including in H d the effect of doublons, through Eq. ((3|), 
within a second-order perturbation expansion, introduces 
[l6j corrections to the transfer integrals, and antiferro- 
magnetic (AF) exchange integrals: 



Ai(R, R') 



t(R)t(R') 

U 



J(R) 



t(R)t(-R) 
U 



(7) 



for R ^ and R' ^ 0. Considering such integrals only 
for nearest-neighbor R and R' results in their expression 
in terms of the parameter J = t 2 /U: 

At(±ax,±ax) = At(±ay,±ay) = At(±ax,±ay) = — J, 
J(±ax) = J{±ay) = J. (8) 

The parameters in Eqs. (I2I8[) have been determined on 
the basis of first-principles calculations [28 30, 36, 37]. 
The set of their values chosen here is: t = 0.43 cV, t' — 
-0.07 eV, t" = 0.03 eV, and J = 0.1 eV (other sets 
of values, within the range of the calculated parameters, 
yield similar results to those below). 

Applying the above approximations, and the con- 
straint in Eq. ([5]) which is approached through a grand- 
canonical scheme, results in the expression of H d in 
Eq. ([1]) in terms of H a , given by fig : 



E 



At(Rj — Rfc, R k — Rj)sl s k ] 



*4-^-.hW}+AH, (9) 

where the mean value A of the constraint Lagrange mul- 
tipliers Xi is treated similarly to the chemical potential 
fi, while A; — A fluctuations are treated through: 



A« = y;(Ai-A)[gk-E a Ua] 



(10) 



H a describes an auxiliary system of fcrmions and 
bosons, behaving according to the laws of physics, but it 
is relevant for the real physical system only under the set 
of values of time-dependent A^ for which the constraint 
is maintained. Assuming lattice translational symmetry, 
let A(k), q^(k) and s£(k) be, respectively, the Fourier 
transforms of Aj, q\ and s\ a to the wave vector (k) rep- 
resentation, determined within the lattice BZ. 



By Eq. (|10p . nonzero Xi — A introduce hybridization 
between QE and between svivon states, and the time 
derivatives of A* induce transitions between such states. 
Consequently, the Lagrangian of the auxiliary system de- 
pends both on A(k), playing the role of generalized coor- 
dinates, and on their time derivatives A(k) which must 
not vanish (for k ^ 0) in order to prevent constraint- 
violating fluctuations. This dependence could be applied 
[HI to derive the conjugate momenta p A (k) of A(k), and 
the dependence of H a on them, reflecting the effect of 
A(k). 

Since A(k) and p A (k) are coordinates and momenta in 
a physical-like quantum system, there exist coefficients 
7(k) for which one can express (for k ^ 0), in a similar 
manner as in lattice dynamics [39j (units where h = 1 
are used): 



A(k) = 7 (k)[J(k) + ft(-k)], 



p x (k) 



2fy(k) 



[/(-k)-/t( k )], 



(11) 



where /^(k) are the creation operators of boson lagron 
states within the BZ (omitting the k = point which 
corresponds to A in Eq. ((9])). Thus, one could express 
AH in Eq. (HHJ) as: 



ah = E 7(q)[Kq) + ' f (-q)] 



(12) 



k,q#0 



x[gt(k)g(k + q)-£4(kK(k + q)]. 

a 

The coupling of QEs and svivons to lagrons, through 
AH in Eq. (fT2")l . introduces an implicit dependence of the 
auxiliary system on A(k) and p x (k), due to the modifi- 
cations introduced to the QE and svivon spectra. The 
effect of A(k) is manifested in the real parts of their in- 
duced self-energies [H, HTJ], and that of p x (k) in the cor- 
responding imaginary parts. These effects on the QE and 
svivon spectra introduce modifications in the Hclmholtz 
free energy T a of the auxiliary system which could be ex- 
pressed in terms of an effective lagron Hamiltonian H x , 
depending on A(k) and p x (k). 

The effect of the magnitudes and rates of the A^ — A 
fluctuations on T a could be generally treated within 
a second-order perturbation expansion. Since a con- 
strained minimum of J- a cannot be lower than the un- 
constrained one, the contribution of each of these time- 
dependent Lagrange multipliers to H x must not be neg- 
ative. Thus, in analogy to the harmonic approximation 
in lattice dynamics [3l|, [3t|, the dependence of H x on 
A(k) and p x (k) has a positive definite quadratic form. 
Consequently, by applying the canonical transformation: 

A(k) = cos(^)A(k)+sin(^)p A (k), 
p x (k) = cos(£ A )/(k)-sin(^)A(k), (13) 

an appropriate choice of results in an expression of the 
form: 



n x s 1 ^[A A (k)At(k)A(k) + J B A (k)p A t(k)p A (k)] 



(14) 
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where the coefficients A A (k) and B x (k) are positive. Ex- 
pressing A(k) and p A (k) through Eqs. (|llll3p yields: 

H x = ^w A (k)[^(k)Z(k) + ±], where (15) 

W A (k) = {[sin 2 (^)v4 A (k)+cos 2 (^)S A (k)] 

x[cos 2 (^)^ A (k)+sin 2 (^)B A (k)]}^ (16) 

lf sin 2 (^)A A (k)+cos 2 (^)i3 A (k) ^ 
171 jl " 2\cos 2 (£ A )A A (k)+sin 2 (£ A )5 A (k) J ' 

The derivation of T~L X from TL a is conceptually anal- 
ogous to the derivation of a phonon Hamiltonian from 
the Hamiltonian of electrons and nuclei in a crystal. In 
the same manner that the auxiliary system represents 
the real physical system only under the correct lagron 
spectrum and coupling constants, a system of electrons 
coupled to phonons in a crystal precisely represents the 
physical system only under the correct phonon spectrum 
and coupling constants. 

The evaluation of the lagron energies cj A (k) in Eq. (fT5")) 
could be carried out self-consistently, on the basis of *H a 
in Eq. ©, applying diagrammatic techniques [39|, H3| 
(this introduces lagron linewidths, ignored in Eq. (|15[) ). 
The coefficients 7(k), appearing in Eq. (TTT1) . should be 
chosen such that Eq. (fT6|) is satisfied, and the svivon and 
QE spectra, obtained on the basis of Eq. (j9|), maintain 
the constraint in Eq. (|6"T). 

These conditions are approached by checking trial 
u; A (k) spectra, determined on the basis of physical 
considerations (see below), and looking for the coeffi- 
cients 7(k) for which these spectra are obtained self- 
consistently. Convergence is reached when the lagron 
spectrum and 7(k) satisfy Eq. (fl"6|) . and determine QE 
and svivon spectra which maintain the constraint; it is 
checked through specific relations between these spec- 
tra, expressed through the "constraint susceptibility" [l6| 
(see below) . It turns out that the derivation of the ma- 
jor features of the lagron, QE and svivon spectra could 
be performed without an elaborate point by point self- 
consistent calculation. 



III. GREEN'S FUNCTIONS 

The elements of the "normal" svivon and QE Green's- 
function [H, Ho| matrices Q_ s and Q_ q are based, in the 
site representation, on expectation values of the form 
( s ie s l<y) an d (<7;9j}j respectively. The Bose condensation 
of svivons, and the pairing of QEs introduce "anomalous" 
[40l | svivon and QE Green's-function matrices T s and 
based, in the site representation, on expectation values 
of the form (s^s^) and {%<lj), respectively. 

The above Green's-function matrices are presented di- 
agrammatically as the propagators shown in Fig. HJa); 
in the k and Matsubara representations [3^, H(| they are 
diagonal and expressed, at temperature T, as functions 
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FIG. 1: (a) Propagator diagrams of the normal and anoma- 
lous svivon and QE Green's-function matrices Q_ 3 , Q_ q , T s , Q_ q , 
and of the renormalized transfer matrix t (see discussion in 
the text); (b) bubble diagrams of the zeroth-order normal and 
anomalous electron Green's-function matrices Qt and £5; (c) 
diagrams of typical first-order corrections to and T^; (d) 
diagrammatical presentation of the geometrical-series sum de- 
termining the self-energy correction E d to G% for [4l[ T > T*. 



of us v = VTtk B T, where fc B is the Boltzmann constant, 
and v is an integer which is even for bosons and odd 
for fermions. The dependence of the Green's functions 
on zcj„ is analytically continued [1^, H(| to the complex 
z plane, including the ±i0 + vicinity of the real w-axis 
(where + is an infmitcsimally small positive number). 

The normal and anomalous electron Green's-function 
matrices Q_ d and T_ are also diagonal in the k represen- 
tation. By Eq. ([5]), the elements of Q_ d (far below the 
upper Hubbard band) are based, in the site represen- 
tation, on expectation values of the form s io .st gt), 

and of F d on expectation values of the form (q i s^q^s^}. 
Thus the zeroth-order (in t/U) normal electron Green's- 
function matrices Q_ Q correspond, in the site representa- 
tion, to bubble diagrams formed by Q_ s and G_ q , and the 
anomalous matrices to bubble diagrams formed by T s 
and as is presented diagrammatically in Fig.fljb). In 
the k and Matsubara representations, and T_o become 
convolutions of such bubble diagrams (see below). 

The relevant term for single-electron excitations in H a 
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in Eq. © could be expressed (using Eqs. <\5\6\i ) as: 

K = ^h s \A s jv<lji where 

ija 

h = ^-(e d + A-M) + (l-^)[t(Ri-Rj) 

+ <jfe Ai ( R « - R k, R k - R,-)L (17) 

where n|. fc is the conditioned expectation value 

(sj. _s fe _ ) when there is a <r electron in site i and no 
electron in site j, or vice versa. 

The elements t.y in Eq. (|17[) introduce a matrix i which 
is diagonal in the k representation, and its eigenval- 
ues correspond (considerably below the upper Hubbard 
band) to an LDA-like band structure, within the one- 
orbital model. This matrix is also presented diagram- 
matically in Fig. [Ija) . The dominant corrections to Q_ d 
and T d , beyond and Tq, are introduced through ?{° 
in Eq. (j 1 T[) ; first-order corrections of this type are pre- 
sented, diagrammatically, in Fig. HJc). 

At T above the QE pairing temperature [41| one has 

F q = Z d = (though F s ^ 0). Similarly to the sum- 
mation of bubble diagrams within the RPA [ljl|4(|, the 
self-energy correction to (due to multiple QE-svivon 
scattering) could then be evaluated, through Eq. (fT7|) . as 
a sum of a geometrical series, presented diagrammatically 
in Fig. [D(d), yielding: 

oo 

E d = ij2(9&r = fa - g$ -1 . (is) 

n=0 

Using Dyson's equation [Uli^, (Q 4 )- 1 = (g£)- 1 -'E d , 
one can express Q_ d as: 

G d = (1 - g$) (1 - (19) 

The consequence of Eq. p^|) is that there are two types 
of poles in the z-dependent Q_ d . One type, contributed 
by the Q c ^ term, consists of a non-FL distribution of con- 
voluted QE-svivon poles; the other type is of poles con- 
tributed by the multiple-scattering (1_ — 2Q^ j t)^ 1 term, 
and it is consistent with electron poles within FL theory. 

Since the evaluation of physical properties can be ex- 
pressed in terms of the electron Green's functions, and 
their poles, these properties are characterized by FL-like 
and/or non-FL features due to the different types of poles 
in Q_ d . The phase diagram of the cuprates [42[ indicates 
a non-FL to FL crossover in the normal-state behavior 
of the cuprates between the underdoped and the over- 
doped regimes 0, [HI- l n heavily overdoped cuprates, a 
crossover is expected to a regime where the large-JJ ap- 
proach, applied here, stops being valid, and there is no 
SC state. 



IV. SPECTRAL FUNCTIONS 
A. Lagrons 

Low-energy (soft) lagron modes play a significant role 
in the lagron spectrum. In the phonon case, the emer- 
gence of a soft mode requires (when it is not at the 
close vicinity of k = 0) a considerable self-energy cor- 
rection, through electron-phonon coupling, in order to 
lower the "bare" phonon energy; by the Kramers-Kronig 
relation (39L l40l | between the imaginary and real parts of 
the phonon self-energy, such an energy lowering means 
that the soft phonon mode has a considerable linewidth. 

On the other hand, in the lagron case, the entire mode 
energy w A (q) consists of a self-energy correction, due to 
its coupling to QEs and svivons, through AH in Eq. (|12[) . 
Thus, the Kramers-Kronig relation of the lagron self- 
energy implies that both the energy and the linewidth of 
a lagron mode could be extremely close to zero. Apply- 
ing this relation, and Eq. (fl4| . indicates that A x (q) and 
-B A (q) almost vanish at this point, q = Q m , of the lagron 
energy minimum, and have the same type of dependence 
on q — Q m within a BZ range of q ~ Q m . 

Consequently, by Eq. (TIB)) . w A (q) also has this type of 
dependence on q— Q m , within this BZ range, and |7(q)| 2 
is fairly independent of q there. This behavior of 7(q) 
is very different from the k dependence of the coupling 
constants between electrons and acoustic phonons in met- 
als which vanish [3l|, for k — > 0. As will be discussed 
further below, the q dependence of w A (q) which is consis- 
tent with the physics of the cuprates is the one yielding a 
Bose-condensation-type behavior. For their 2D BZ, this 
corresponds to w A (q) oc |q — Q m |, for q ~ Q m ; the ap- 
parent singularity of w A (q) at q = Q m is a self-consistent 
consequence of the macroscopic value of (?'(Q m )Z(Q m )}. 

By Eq. (TlTj) , the Q m points of the lagron spectrum cor- 
respond to particularly large values of A(q) and p A (q), for 
q = Q m . This is typical of wave vectors corresponding 
to major spin and charge fluctuations, where an inho- 
mogeneous enforcement of the auxiliary-particles' con- 
straint is essential. Theoretical studies [43l - l45j . based 
on Hamiltonians like 'H d in Eq. {T]), predict that the in- 
terplay between the effects of electron hopping and AF 
exchange is likely to drive the formation of stripe-like in- 
homogeneities (studied in previous works by the author 
[H, Hll ) . Thus the physics of the cuprates is consistent 
with Q m points which correspond to the wave vectors of 
the striped structures. 

In agreement with the above theoretical predictions, 
the hole-doped cuprates are characterized by dynamical 
or static stripe- like inhomogeneities [l8l-|22|: in SC sto- 
ichiometrics, they correspond to stripes directed along 
the a- or the 6-axis; however, an intrinsic mechanism (of 
a nature discussed below) introduces long-range symme- 
try between stripe segments directed along these axes, 
resulting in a checkerboard-like structure [22[ ■ Con- 
sequently, the lagron spectrum corresponding to such a 
structure is of the type presented in Fig. [2J 
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FIG. 2: A typical lagron spectrum in hole-doped cuprates; the 
minima correspond to degenerate striped structures p^j - [2^ |. 



merit of the analogous electron-phonon coupling term in 
phonon-mediated SCs is based on the fact that the ratio 
between phonon and electron e nerg ies is very small; con- 
sequently, by Migdal's theorem [25|], "vertex corrections" 
beyond a second-order diagrammatic expansion are neg- 
ligible due to a small integration phase space. 

Concerning the treatment of AH, as can be viewed 
in Fig. [2 the lagron spectrum consists of two regimes; a 
low-energy one around point Q (including the points Q m 
defined in Eq. (|2T))) ). and a high-energy regime in the rest 
of the BZ. As is discussed below, the peculiar behavior 
(including pairing) introduced by AH is primarily due 
to lagrons of the low-energy regime, for which Migdal's 
theorem applies. Vertex corrections due to lagrons of the 
high-energy regime may be not negligible, but they are 
regular and could be approximated as renormalization 
effects within the second-order diagrammatic expansion. 



Thus, the lagron energy band, w A (q), is characterized 
by V-shape minima: 

uj x (Q m ) = k B T/0(N), at 

Qm = Q + Sq m , for m = 1,2,3 or 4, (20) 

where Q = (% / 'a)(x + y) is the wave vector of the AF or- 
der in the parent compounds; <5q m = ±Sqx or ± Sqy are 
modulations around Q, indicative of the striped struc- 
tures [H|-[22j (typically, Sq = ir/Aa). The lagron spec- 
trum presented in Fig.[2]corresponds to degenerate states 
of stripes directed along the a and b axes, and of differ- 
ent phases of the spin periodicity (see below). The four 
lagron energy minima, and the symmetry of the system, 
determine the main features of the spectrum, and as was 
discussed above, the linewidths are generally larger at 
higher energies. The high-energy extent of this spectrum 
is somewhat higher than that of the svivon spectrum [l6| , 
and is thus expected to be ~ 0.3-0.4 eV. 

The existence of four equivalent lagron minima in 
points Q m in Eq. (|20[) results in a combination state of 
four Bose condensates of the lagron field [46| , each con- 
sisting of a striped structure of the symmetry correspond- 
ing to one of the Q m wave vectors. In SC stoichiome- 
trics, such a combination state is energetically favorable, 
for T — > 0, on one of the symmetry- broken condensates 
(corresponding to a long-range striped structure) since it 
is the one which yields pairing (see below). Similarly to 
phonons which are Goldstone bosons [47j associated with 
spontaneous breakdown of lattice translational symme- 
try, the lagrons are Goldstone bosons associated with the 
breakdown of the symmetries of the long-range striped 
structures, corresponding to each of the combined con- 
densates. 

The svivon and QE spectra are evaluated through 
a self-consistent second-order diagrammatic expansion, 
where a mean-field treatment of H a in Eq. ([9]) is applied 
at the zeroth order. The expansion is carried out on two 
Hamiltonian terms; one of them is the svivon-lagron and 
QE-lagron coupling term AH in Eq. (Ti"2"|) . The treat- 



B. Svivons 

In the parent compounds, the lagron spectrum includes 
one V-shape minimum at q = Q, corresponding to a 
condensate (see above), and the Bose condensation of 
the svivon field [46| is manifested in long-range AF order 
[33l | . Within such a svivon condensate which corresponds 
to a combination of Neel states, shifted from each other 
by ax, the expectation value n s = {sl a s ia ) is site- and 
spin- independent (and equals 0.5 for zero doping). 

The svivon band is then characterized [33| by V-shape 
minima of energy k B T/0(N) at the points ±k = ±Q/2. 
There are four inequivalent values of ko, introducing a 
broken symmetry: 

k, S §=±£(* + tf) «*£(*-<>). (21) 

They correspond to two degenerate condensates of rect- 
angular symmetry, with axes along the x + y and x — y 
directions. Within the grand-canonical scheme, the 
svivon minima at ±ko in the AF condensates are self- 
consistently formed due to [33[ large values of both 

n s (k) = (4(kK(k)}, and m s (k) ee (^(k^-k)) 
(note that n s = — ^ n s (k)), (22) 

k 

around the minima (becoming O(N) at k = ±ko). 

In SC stoichiometries, the minimum of the lagron spec- 
trum at point Q splits into the four Q m points, intro- 
duced in Eq. (|2"0|) . and shown in Fig. [5] As was men- 
tioned above, this spectrum corresponds to a combina- 
tion of four condensates, in which AF order is replaced by 
striped structures, corresponding to svivon condensates 
[46| of energy minima which are shifted [l6| from ±ko 
(see discussion below) . 

The svivon minimal energies, within their condensates, 
are then [16[ not as close to zero as k B T/0(N) (and the 
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shape of the minima becomes parabolic close to the bot- 
tom), since they are evaluated on the basis of a lower- 
free-energy state (due to coupling to the QEs) which is 
a combination of the condensates (and thus it consists of 
fluctuating striped structures); but one still has (within 
each condensate) m s (k) ^ which are large, together 
with n s (k), around the energy minima (though not as 
large as 0(N) — see below). 

Note that A(q) and p A (q) in Eq. (fTTj) scale with 
Z(±q) and ^(±q) which behave like (^(±q)/(±q)) 2 for 
q — Qm , while n s (k) and m s (k) are determined through 
Eq. (HU)). Thus, the fact that the svivon energy minima 
do not get too close to zero enables such a lagron spec- 
trum (thus with the V-shape energy minima specified in 
Eq. (HO))) to yield values of A(q) and p A (q) which have 
appropriate values to maintain the constraint 46] within 
the q ~ Q m ranges. 

Bose-condensed systems are characterized by nonzero 
expectation values of the field operators [39|, [4(| . Treat- 
ing the (combined) svivon condensates adiabatically, one 
can express, through Eq. (j2"2")l . m s (k) = (s^.(k))(s, (— k)); 

these (s*(k)) and (s,(— k)) have the same magnitude, 
but different phases due to the different sites where the 
contribution to a =j" and a =], svivons is maximal, and 
their phases fluctuate in a correlated manner. Terms in- 
cluding such expectation values in the final result of a 
calculation, do not vanish due to phase fluctuations only 
if they appear as products of the form (s^(k)) (s, (— k)) 
(which are expressed as m s (k)). 

The second Hamiltonian term, treated through a 
self-consistent second-order diagrammatic expansion, ac- 
counts for the coupling between QEs and svivon fluctu- 
ations which are expressed through terms of the form 
s CT (k) — (s CT (k)). It is obtained by expanding H a , in the 
rhs of Eq. ((9|) , in such fluctuations (when the svivon op- 
erators there are expressed in the k representation), and 
treating them as a perturbation. This Hamiltonian term 
could be approximated as (see Eq. (IT71) '): 

n ' ~ ^7 E {?(k' + k")(4(k' + k"-k)) g t( k ) 

kk'k"<x 

xg(k')[4(k")-(^(k"))]+^.}, (23) 



where t(k) is the Fourier transform of t(R) which is, ap- 
proximately, defined similarly to t(R) in Eq. ([2]), mod- 
ifying t' and t" (due to the At terms in Eq. ((8|)) to 
£' = £' + 2n s J and t" = t" + n s J, respectively (while 
t remains unchanged). 

Following the standard treatment of Bose-condensed 
systems [33, H(| , the (Bogoliubov) boson creation opera- 
tor of an excited svivon state of wave vector k and spin 
a is derived as [33j a combination of sj.(k) and s_ a (— k). 
A 2 x 2 energy plus self-energy matrix is derived within 
the basis created by the spinor operator (s|(k), s^(— k)). 
It is expressed as 

eg(k) + £ s (k,z) = [e s (k)+^ s (k,z)}l + m s (k,z)r 1 
+3$ s (k,z)r 2 , (24) 
where 1, t_ y and r 2 are the unity and Pauli matrices, and 
z are "complex energies" [33, |40|. The diagonal terms 
in the rhs of Eq. (|24p include £g(k), derived within the 
mean-field approximation, and S s (k, z), derived through 
the diagrammatic expansion. The non-diagonal terms, 
$ s (k, z) and $ s (k, z)* , are derived through both the 
mean-field and the diagrammatic calculations; they exist 
due to the nonzero values of m s (k'). 

Two degenerate Bogoliubov svivon states, correspond- 
ing to (k f) and (— k 1), are derived as a function of 
z, and their boson annihilation operators are denoted by 
s^.(k, z) and s^(— k, z), respectively; s^(k) and s^(— k) 
are then expressed as [33$ : 

S >k) = exp(^( Z )/2)[cosh(^(z))5 CT (ak,0) 

+ sinh(^(z))sL CT (-ak,z)], (25) 

where the coefficients ££(z) and </£(z) are determined 
by the requirement that the matrix eg(k) + S s (k, z) in 
Eq. (|2"4")l is diagonalized under the transformation to 
these Bogoliubov states. This transformation diagonal- 
izes Dyson's equation (33") Eq| , yielding the poles of the 
diagonalized svivon Green's function Q s (k,z) at: 



z = e s (k) + £ s (k,z), where Q=£ s (k, z) = cosh (2££(z))9£ s (k, z), 



c5(k) + 5ffi s (k, z) = E s (k, z) = sign[eg(k) + 5RS s (k, 2f)] J[eg(k) + 5R£ s (k, z)] 2 - | cos (^(z) - </£(z))$ s (k 



$ s (k,z) = |$ s (k, z)\exp(iipk{ z ))> an d | cos (ip^{z) - <fc(z))\ = min 

I 



1, 



eg(k) + KE s (k, z) 



$ s (k,z) 



(26) 



The coefficients ££(z) are obtained through: 
eg(k) + 5RS s (k,z) 



The svivon spectral functions are obtained through: 



cosh(2£(z)) 
sinh(2^(z)) 



cos(^(z)-<^(z))|<l> s (k,z) 
E s (Kz) 



(27) 



A s (k,io) ee %g s (k,io -iO+)/n 

r s (k,w)/27r 



[u - E s {k,u)} 2 + [ir s (k,w)] 



r, (28) 
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where r s (k,w) = 23E s (k, w-iO + ). G s (k, z) has positive- 
and negative-energy poles [l6| (see below); r s (k, w) and 
A s (k,Lu) are positive for oj > and negative for cj < 0, 
and the w-integral of A s (including the contributions of 
the different poles) is normalized to one. 

The values of n s (k) and m s (k) in Eq. (|22|) are obtained 
through [H: 

n s (k) = y dc. cosh (2^H)A s (k, W )[& T H + 

m s (k) = J duj sink {2C k {uj))A s (k,uj)[b T {u J ) + |] 

xexp(i#(w)), (29) 

where & T (w) = l/[exp (co/k B T) — 1] is the Bose distri- 
bution function; one has 6 T (w) + | = isign(cj) in the 
high-|o;|/fe B T limit. 

A lagron spectrum of the type presented in Fig. [^deter- 
mines eight degenerate svivon condensates with energy 
minima at the points: 

±^ = ±(k + ^), (30) 

for each of the four values of m in Eq. (|20l) , and the two 
possibilities for ±k in Eq. (|2"Tj) . Lagrons at the ±Q m 
points (of the energies in Eq. (|2"0")0 induce symmetric 
positive- and (lower weight) negative-energy svivon spec- 
tral branches [lfj, corresponding to the splitting of the 
Green's functions poles, due to the inhomogeneities. The 
phases V , k( w ) an d ^kC^) tnen satisfy: 

^-kH = ^k(-w) = + 7T, 

0- k M - = (3i) 

Thus, by Eqs. (|26l3ip . the svivon condensation order pa- 
rameter <£> s (k, uj) reverses its sign upon the sign reversal 
of cither k orw. 

The effects of entropy at high T, and of pairing at low 
T (see below), drive the system to a combination state 
of the eight broken-symmetry condensates, specified in 
Eq. (|30|) . at high T, and of four of them, corresponding 
to one of the two possibilities for ±ko in Eq. (|2"Tj) . at 
low T. Such combinations reflect fluctuations between 
the condensates. Consequently, the planar symmetry of 
the system remains square at high T, while at low T it 
is broken to a rectangular symmetry with axes along the 
x+y and x — y directions, in agreement with experiment. 
Svivon spectral functions in the SC state have been pre- 
sented elsewhere fig. 



C. Unpaired QEs and electrons 

1. General features 

The QE spectrum is evaluated treating fluctuations 
between the combined svivon condensates adiabatically. 
Thus, in analogy to the effect of lattice dynamics on 



single-electron states in a crystal, the contribution of 
coupling to the different svivon condensates to the QE 
states is averaged out, and Q_ q remains diagonal in the k 
representation (corresponding to the lattice plus a super- 
structure determined by Q — see below), though scatter- 
ing is introduced. This conclusion remains approximately 
correct also under the random long-range distribution of 
static inhomogeneities, observed in a certain range of the 
phase diagram [2l|, [22[ , since the variation that they in- 
troduce to the charge distribution, and thus to the effec- 
tive potential of the QEs, is rather small and insufficient 
to introduce Anderson localization. 

Since, by Eqs. (|21l30p . the points ±Q m /2 of the svivon 
energy minima within their condensates [16J, are sym- 
metrical around ±ko = ±Q/2, the averaged svivon spec- 
tral contribution to the QE states is symmetrical around 
these points as well; consequently, it is invariant [l6| un- 
der a k-shift of Q (note that 2Q is reciprocal-lattice 
vector), and thus it corresponds to an effective AF or- 
der. The symmetry of the electron spectrum is the 
same as that of the convoluted QE-svivon spectrum (see 
Eqs. (|5ll9p V Thus, since the electron spectrum is invari- 
ant under a k-shift of 2Q, the QE spectrum must also be 
invariant under a k-shift of Q. 

By introducing appropriate phase factors to the QE 
and svivon states, their k values can be shifted by ko or 
— ko, so that corresponding QE and electron bands ap- 
pear in the same BZ areas. Thus, by Eq. (f30|). a shift of 
— ko in the svivon k values results in the shift of the en- 
ergy minima in their condensates from ±Q m /2 to <5q m /2 
and Q — <5q m /2. 

The evaluated QE spectral functions A q (k, ui) for hole- 
doped cuprates, corresponding to n s — 0.42 (thus close 
to "optimal stoichiometry" ) , at fc B T = 0.01 eV, within 
the "hump phase" [Uj], are presented in Fig. (3Ja-c). The 
calculation is based on a typical svivon spectrum, where 
the averaged svivon-condensate spectra [1 61 ] have broad 
energy minima which are approximately linear in T in 
this phase (see below). 

The determination of the svivon and QE linewidths is 
discussed below, and they are, approximately, consistent 
(within the hump phase) with the marginal-Fermi-liquid 
(MFL) phenomenology [13]. The fact that the QE spec- 
trum is invariant under a k-shift of Q, is reflected in the 
existence of equivalent "main" and "shadow" QE bands 
obtained from each other by a Q shift. 

As was discussed above, the QE spectrum has a mean- 
field "bare-band" part, eg(k), an d a "dressed" part, in- 
troduced by self-energy corrections E 9 = Sf + £?, due 
to coupling to svivon fluctuations (through %') and to 
lagrons (through AH), respectively. The QE spectral 
functions are obtained through: 

A 9 (k,u;) = %g q {k,Li-i0 + )/TT (32) 

r9(k,a;)/27r 
= [oj - e«(k) - 5RS?(k,w)] 2 + [r?(k,w)/2] 2 ' 

where r«(k,w) = 23S]' 3 (k, w-i0+). The effect of inhomo- 
geneity is that the QE Green's function (? 9 (k, z) has more 
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than one pole (see below). The w-integral of A q (k,uj), 
including the contributions of the different poles, is nor- 
malized to one. 

The effect of the nearest-neighbor transfer integral t 
on Cg(k) drops out (due to opposite-sign contributions of 
svivon states shifted by Q). The effects of t' and t" (see 
Eq. P]>) on £g(k) are renormalized, due to svivon and At 
terms in Eqs. (I8l9p . and they could be, approximately, 
replaced by F = 2n s (t'+4n s J) and t" = 2n s (t" + 2n s J), 
respectively. The n s factors here are derived from the 
summation of expectation values of products of svivon 
operators which could be expressed both in terms of the 
?i s (k) and of the m s (k) factors (which have a major con- 
tribution at the vicinity of the svivon energy minima), 
and applying Eqs. {22]), {27]) and ([29]) . 

The independence of e^k) on t reflects the fact that 
a QE represents [l6| an approximate electron for which 
the creation operator of its svivon component, through 
Eq. {5j, is replaced by its expectation value. Since av- 
eraging the svivon spectrum over its degenerate conden- 
sates corresponds to an effective AF order (see above), 
the dependence of the QE spectrum on nearest-neighbor 
hopping processes (which disturb this order) must involve 
svivon fluctuations, through TL' in Eq. {23]) . Thus, the 
effect of t on it is introduced by the self-energy term: 



^( k >-) = ^E 

k'k" 



doj q A q (k',uj q ){ \t(k' + k" 



x |m 5 '(k' + k 
x {cosh 2 (^(^)) 



k)| j J du s A s {k" ,lo s ) 
b T (u J s ) + / T (-w«) 



+ sinh 2 (^( W s )) 
|m s (k" 



z — lo s — ojQ 
b T (u s ) + f T (oji) 



z + ui s — (jjl 
,"M2| 



]} 



|i(k + k")| 2 |m s (k + k"-k')| 

z — <jj1 J 

dwM s (k",w s ){cosh 2 (^(a; s )) (33) 

+ sinh 2 (^(a/)) 
|m s (k"^ 



b T {u s ) + f T (oj") 



z + lo s — ojQ 

b T (LU S ) + / T (-0 




Z — U! s — U) q 

where f T (co) = l/[exp (w/fc B T) + f] is the Fermi distri- 
bution function. 

By Eq. {2"9"]l . the magnitudes of the |m s (k")| terms in 
the rhs of Eq. ([33]) . and of the integrated terms that 
they are subtracted from, differ (in the vicinity of the 
svivon energy minima, where the major contribution to 
the k" summation comes from) mainly because of their 
different energy denominators. When the signs of the 
denominators of these terms are opposite, they both con- 
tribute to "pushing" the QE energies towards zero. Con- 
sequently, the 5^-induced self-consistent renormalization 
of the QE energies is a shift towards zero which becomes 
larger closer to zero. 



This results in BZ areas of "flat" QE bands which 
would have been at energies of the k B T scale, if it were 
not for the effect of QE-lagron coupling (introducing to 
them the structure discussed below). The transition be- 
tween these low-energy areas and the higher-energy BZ 
areas is through an almost discontinuous wide-energy 
zone, as is viewed in Fig. (3Ja-b). This peculiar spectral 
structure reflects the distinction between strongly renor- 
malized low-energy QEs which are subject to hopping 
processes maintaining the stripe-like inhomogeneities, 
and high-energy QEs which are not; it is primarily deter- 
mined by the contribution of the hopping parameter t to 
W in Eq. {23]). 

The corresponding electron spectral functions 
A d (k,uj) = ^sQ d (k, io — iO + )/ir, obtained on the basis of 
Eq. {T9]) . are presented in Fig. EJd-f). The high weight 
of the svivon poles around their minimal energies [161 ] 
results in band-like features in the electron spectrum 
which follow the QE band (see below). 

The effect of the multiple-scattering FL-like electron 
poles in Eq. (|19p is the removal of the equivalence be- 
tween the main and the shadow electron bands, as is ob- 
served in experiment (48, 49]. However, non-FL behavior 
associated with the QE Fermi surface (FS) in Fig. (3Jc) 
persists, contributing the hole and electron pocket-like 
features appearing in Fig. [3jf), and detected in experi- 
mental data (5p| mainly in the underdoped regime. The 
existence of such FS pockets has also been confirmed in 
the observation of quantum oscillations (5^, [Hj] . On the 
other hand, the FL-like contribution to the normal-state 
electron spectrum is dominant in the overdoped regime, 
where quantum oscillations typical of a large FS have 
been observed [53] ]. 

The low-energy flat QE bands, and the almost discon- 
tinuous transition from them to the high-QE-energy BZ 
areas are modified in the electron spectrum (see Fig.[3Jd- 
e)) to low-energy kinks followed by "waterfalls" , in agree- 
ment with experiment (5~il. [55|. As was discussed above, 
this spectral anomaly results from the renormalization of 
the QEs, such that their low-energy excitations maintain 
the stripe-like inhomogeneities. 



2. Low- energy range 

The self-energy term, introduced by QE-lagron cou- 
pling through AH, in Eq. (|f 2[) is expressed (for un-paired 
QEs) as: 

£*(k,z) = ^E/ d^A^k-^uj") 



F T (q,u*,z) 



xF T (q,uj q ,z), where (34) 
"M^(q)) + / T (-^) 



I7iqj 



z — 0J x (q) — oj q 
b T (u x (q)) + f T (cj q ) 



z + uj x (q) — ui q 
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FIG. 3: The QE (a-c) and electron (d-f) spectral functions in optimally hole-doped cuprates along a "line of nodes" (a,d), along 
a parallel line (b,e) shifted by 0.4257r(y — x), and on the "Fermi surface" (c,f), determined at u — — 0.1fc B T. 



Within the low-QE-energy BZ areas, the w-dependence 
of is weak, and its imaginary part is relatively small. 
Consequently eg could be renormalized there to include 
the effect of Ef , and the anomalous behavior of A q {k, lu)) 
is obtained approximating E 9 in Eq. ([32]) by the QE- 
lagron term E* in Eq. 

As is viewed in Fig. [21 one could specify the lagrons 
according to five ranges in the BZ: (1) those at the four 
energy minima points Q m in Eq. (|20p , referred to as "Q m 
lagrons" ; (2) those at the lagron "extended saddle point" 
(ESP) around point Q, referred to as "Q-ESP lagrons" 
which introduce an energy scale ~cj a (Q); (3) those of 
energies < w A (Q), around the Q m points, referred to as 
"Qm-vicinity lagrons" ; (4) those around the high-energy 
saddle points at (n/a)x and (n/a)y, referred to as "high- 
energy SP lagrons"; (5) those within the rest of the la- 
gron spectrum, referred to as "continuum lagrons" which 
occupy most of their phase space. 

Analogous expressions to Eq. (IMl) . where A q , u q and 
f T (±uj q ) arc replaced by A s , lu s and — b T (±uj s ) (and 
cosh 2 and sinh 2 factors are introduced — see below), are 
obtained for the contribution of AH in Eq. (fT2")) to the 
svivon self-energy terms, appearing in Eq. ([24[) . The 
Q-ESP lagrons then help stabilize high-spectral-weight 
svivon states somewhat below w A (Q) in the SC state. 
Consequently jlc| . the resonance-mode energy E res is 
somewhat above cj a (Q) [56]. 

The spectral functions of low-energy QEs, coupled to 
each other through Q TO lagrons, split (due to the large 
values of the Bose functions in Eq. ([34)) — see Eq. ([20[) and 
the discussion above) into positive- and negative-energy 
peaks, referred to as "humpons" ; such a splitting is ex- 
pected due to the spin-density waves (SDW) associated 



with the stripe-like inhomogeneities which are generated 
by the lagron spectrum in Fig. [5] 

As was discussed above, within some range of temper- 
atures, free energy (where the effects of both energy and 
entropy are accounted for) is minimized for a phase of 
fluctuating inhomogeneities. For unpaired QEs, a third 
QE spectral peak appears in this phase, between the 
positive- and negative-energy humpons (as is sketched 
in Fig. IUa)), and it is referred to as a "stripon"; it rep- 
resents low-energy charge carriers due to the fluctuating 
charged stripes. 

A stripon at point k is associated with the pole of 
Q q (k, z) which corresponds to lattice periodicity, while 
the humpons result from poles corresponding to the in- 
homogeneities. The energies of the humpon peaks 
correspond to the energy associated with the inhomo- 
geneities, and could be as high as ~J. As is discussed be- 
low, the stripon peak is sharper than the humpon peaks 
for S; B TCu H , and its width increases with T. 

The 7(q) coefficients in Eq. (TT21 remain finite for q — > 
Q m (see Eq. ([To) and the above discussion) and their 
values there, as well as the values of w A (Q m ) in Eq. ([2TJ]) . 
are self-consistently determined to yield such humpon en- 
ergies due to coupling to lagrons, and specifically those 
around the Q m points, through Eq. Q34p. The contribu- 
tion of the Q-ESP lagrons results in w H >^w A (Q). 

Within the low-QE-energy BZ areas there are "nodal 
areas" , close to the lines of nodes y — ±i, and "antinodal 
areas", close to circles of radii of about (n / 4a)-(n /3a) 
around the (ir/a)x and (w/a)y points. As can be seen in 
Fig. \3lc), a difference between the QEs in the antinodal 
and nodal areas is that the first, but not the latter (if they 
are close enough to the lines of nodes) , can be coupled to 
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FIG. 4: Typical low-energy QE (a) and electron (b) spec- 
tral functions in the normal and SC states, in nodal (ar- 
con, Fermi arcs) and antinodal (humpon, stripon, peak, dip, 
hump) points in hole-doped cuprates. Their structures are 
presented in the MFL, PG, and SC states (see discussion in 
the text). 



other low-energy QEs through (up to four) Q m lagrons 
(see Eq. pj|)). 

Consequently, the spectral functions of unpaired QEs 
within the nodal areas (referred to as "arcons" ) are char- 
acterized by a simple peak (due to one-Q q pole per 
k point), of the k B T energy scale, as is sketched in 
Fig. Eta). On the other hand, the spectral functions 
of QEs within the antinodal areas are characterized by 
the above humpon-stripon-humpon structure. As will 
be discussed below, these spectral features of the QEs 
(specified as arcons, stripons and humpons) are associ- 
ated with different physical features in the cuprates. 

Thus a two-component scenario is described here, con- 
cerning the symmetries of the low-energy QEs in the 
nodal and antinodal areas, due to the different number 
of poles in {? 9 (k, z) they correspond to. An increase in 
T results in an increase in the spectral weight within the 
stripons peaks in the antinodal areas (on the expense of 
those within the humpon peaks), and in the shift of the 



borderline between the nodal and antinodal areas (con- 
cerning the symmetry of the QEs) towards the antinodal 
ones. When k B T crosses w H , a crossover occurs from the 
above inhomogeneity-derived two-component scenario to 
a homogeneous one, where all the low-energy QEs corre- 
spond to one-Q q pole per k point. 

For k B T <^ ui H , the arcon and stripon peaks at different 
k points can be treated as an almost flat band. The 
energy center of this band is close to the centers of most of 
the peaks, determined through (see Eq. (|32|) ) the relation: 
ui = eg(k) + K£«(k,w). The T dependence of the band 
center approximately corresponds to that of the effective 
QE chemical potential /j, — A (see Eq. ©); thus this 
T dependence scales with k B T, if the band is not half 
filled and its width in not larger than ^k B T, and it is 
weaker, otherwise. 

The above low-energy spectral features of the QEs are 
consistent with the nodal-antinodal dichotomy occurring 
in the cuprates. Furthermore, they provide convenient 
QPs to study, discuss, and analyze (see below) their pe- 
culiar low-energy spectrum and its consequences on the 
derivation of their anomalous physical properties. 



3. QE scattering rates 

The scattering rates T«(k,w) = 23E«(k,w - i0+) of 
unpaired low-energy QEs can be expressed, on the basis 
of Eq. dM|), as: 

r«(k,u,) - ^]T| 7 (q)| 2 {^(k-q, W -,/(q)) 
q 

x[6 r ( W A (q)) + / T (^(q)- W )] 

+A«(k-q,w + w A (q)) 

xMw A (q)) + / T ( W A (q)+ W )]}, (35) 

where an approximately q-independent |7(q)| 2 could be 
assumed (see Eq. and the above discussion). 

Eq. (1331) is applied to study the behavior of T q , specif- 
ically in the low- and high-|w|/A: B T limits. As can be 
viewed in Fig. [3Ja-b), the k q points where A q con- 
tributes to Eq. (1331) . for < 0.3 eV, consist mainly 
of those within the low-QE-energy BZ areas, shown in 
Fig-EJc). The lu dependencies of A q (k, uj) within the low- 
energy arcon and stripon peaks, the mid-energy humpon 
peaks, and the high-energy areas, are specified in terms 
of canonical functions A q c (uj), A^^lo) and A^ c (w), re- 
spectively. 

The energies at the maxima of A^^uj) are generally 
too high to have a significant contribution to T q , through 
Eqs. (|32|35j) . within the ranges of T and ui studied here. 
The humpon peaks (characterized by A^ e ) exist within 
the low-energy antinodal areas in the fc B T « w H regime, 
and if also \uj\ <C w H , their role is negligible, similarly 
to that of A^; for |w| > uj h , they are approached as a 
case of A? . In the k B T > u H regime the humpon and 
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FIG. 5: Typical results (in arbitrary units) for the QE scatter- 
ing rates T q , at [41[ T > Tg, in k points ranging between the 
lines of nodes and the antinodal BZ areas, as (a) a function 
of T for |w| < fc B T p 9 , and (b) a function of w, for fc B T < rg a 
(see discussion in the text). 



of T and w, in the low \cj\/k B T limit, and at T above [4J| 
T« and below T q Jk s (see Appendix A) and w A (Q)/fc B , 
respectively; these results include k points ranging be- 
tween the nodal and antinodal BZ areas. Approximate 
linear dependencies of r 9 (k, lu) on T and on u are found 
in the low- and high-|w|/fc B T limits, respectively, above a 
k-dependent low-energy scale u; L (k). The high-|w extent 
of this linear dependence is limited by the high-energy ex- 
tent of the lagron spectrum in Fig. [2] (which could exceed 
~0.3 eV), and the high-T extent is limited by the phase 
stability. The contribution of humpons to T q results in 
the increase of the slopes of the T q vs k B T and \co\ curves 
when their increasing values approach w H . 

For k B T & \oj\ <C w L (k), r q (k, uS) can be approximated 
1000 as Tg(k) (determined in Eq. (|A3[) ) which is close zero on 
the lines of nodes, and to in the antinodal areas. As 
is discussed in Appendix A, the value of w L (k) is close 

and is somewhat smaller 

HQ) 



to Tq, in the antinodal areas 



than 1 5 



l (Q) 



\E on the lines of nodes; 



is smaller than Fg a in the heavily underdoped regime, 
and their values cross each other when the doping level 
is increased. 



4- Svivon scattering rates 



due 



m 

, ( , due to svivon-QE 
has a low-energy con- 



The svivon scattering rates V s include a term T s x , 
to their coupling to lagrons through AH in Eq. ([1 
yielding an analogous expression to that for T q 
Eq. (|55|) . However, also a term T 
coupling, through %' in Eq. (|23 
tribution to V s which becomes significant above [4l| T q . 
Since A s (k, lj) includes symmetric positive- and negative- 
energy branches (of different weights), and the sign of 
r s (k, (jS) is that of to, one has (differently from QEs — see 
Eq. (fA3]P l: 



stripon peaks are merged into an arcon peak which is 
again approached as a case of A q e . 

The evaluation of r 9 (k,cj) is detailed in Appendix A; 
typical results are presented in Figs. EJa-b), as functions 



T s (k,io -+0) =0. 



(36) 



Using Eqs. (JUJ), (H3J) and (gH), one can express V s (for 
unpaired QEs) as: 



T s (k, W ) = r^(k, W )+r^(k,w), where 

rj(k, W ) S 2 ^ C ° Sh ^ H) E I^(k + k")| 2 |^(k + k"-k0| j du'A%V,J)A%V\J-u)[f T {J-u)-f T ("')] 



k'k" 

S ^!|iMl ^ | T(q) p( #(k _ q _ w _ w A (q)) cosh 2 (S _ q(w _ ^(q)))[6 r (^(q)) - fc T (^(q) - w )] 
q 

+A S (k - q, u + w X (q)) cosh 2 (g_ q (w + c A (q))) [b T ( W A (q)) - b T (lu X (q) + «)] } (37) 
+A*(k - q, c A (q) - w) sinh 2 (^_ q (c A (q) - c))[6 T ( W A (q) - U ) - M^(q))] 

(k - q, -u> - ui X (q)) sinh 2 (££_ q (-c - CJ X (q))) [6 T ( W A (q) + lu) - b T (lu X (q))] } . 

I 

Let e s (k) be the positive svivon energies [l6| (within are the lower- weight negative energies — see above); let 
their condensates) at the maxima of A s (k,Lu) (— e s (k) 
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^min = e s (k m in) be the minimum of e s (k) (within each 
condensate there are two k m j n points, separated from 
each other by a Q m vector). Above [4l[ Tg, the T de- 
pendence of e s (k) is determined by the variation of A in 
Eq. ([9]) to maintain a T-independent n s through Eqs. (f22l 
[29]) . By Eq. (|26|) . this results in the approximate scaling: 



e s (k) oc v^f" 



for Ik — k n 



kjnin | ) 2 4 



(T-T ) 2 , 



(38) 



where c s and To are constants (one has To < and 
it could be in principle either positive or negative — see 
below), and A; ma x represents the radius of an approximate 
circle around k m i n , in the BZ, where there is a noticeable 
effect of <& s on I s (k). 

By Eq. (f38|) . e^ in approximately scales with T — T , 
such that it remains sufficiently larger than k B T (as 
[l6| below T^) that 6 T (|u;|) is relatively small (though 
not negligible), within the range where k ~ k m ; n , and 
|^4 s (k, a;) | is significant. On the other hand, Eq. (f38|) 
yields that the effect of T on e s (k) ^> e^ in is through 
a minor additional term scaling, approximately, with 
(T — T ) 2 , and in such k points one has 6 r (|w|) <C 1 
within the u> range where |A s (k,cj)| has a considerable 
magnitude. 

When T is increased, the minima of e s (k) around 
e* lin become broader and flatter, and since |T /T| is de- 
creased, the ratio e^ in /T changes as is necessary in order 
to maintain a constant n s , through Eqs. ([221 |2"9")) (un- 
der the opposing effects of T on it through the minima 
and cosh (2££(w) — see below — as well as the effect of the 
linewidth of e s (k), for k ~ k m i n ). 

By Eqs. ((HD, J23) and fl3SJ|, the cosh (2££(w)) factors, 
appearing in Eq. (|37l) . can be approximated as: 



cosh(2£M) 



(c s k max ) 2 + (T - To) 2 



y (c s |k kmin 

for Ik - k min 



|) 2 + (T-T ) 2 ' (39) 
< fc max , and \u\ < e s (k). 



These factors are largest in k points at the low-svivon- 
energy (LE) BZ areas, close to the k m j n points, where 
r s (k, oj) is significant in the \oj\ < e^ in range, and by 
Eq. (|39)l they, approximately, scale there with 1/(T — T ), 
for |w| < e s (k). On the other hand, in k points at the 
mid-svivon-energy (ME) BZ areas, where r s (k, u) is sig- 
nificant in the ~ef nin < \uj\ <~2J range, and Eq. (|3"9")l 
is still valid, it yields a weak T dependence for the 
cosh(2££(w)) factors, through a term cx(T - T ) 2 . The 
cosh (2££ (w)) = 1 limit specifies k points at high-svivon- 
energy (HE) BZ areas, where the effect of condensation 
is missing and Eqs. (|38l [39]) are invalid; T s (k, u) is sig- 
nificant there in the high-w and -T limits. At the LE and 
ME areas one has cosh(2££(w)) = 1 for \u\ > e s (k). 

Thus, one could classify the svivon spectrum within 
the BZ according to LE, ME and HE areas, and crossover 
areas which could be, approximately, split between them. 
Since the minima of e s (k) around e^ lin become broader 
and flatter when T is increased, the LE areas are growing 
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FIG. 6: Typical results (in arbitrary units), at [U T > T$, 
for (a) the T dependence of the to derivative rg' of the svivon 
scattering rates F s at oj = 0, in LE and ME svivon k points 
(see discussion in the text) and of its k-integrated value, and 
for (b) the w dependence of T a , in LE, ME and HE svivon k 
points. 



with T on the expense of the ME areas. By Eq. (|38p . the 
size of the LE areas scales with (T — T ) 2 , and thus the 
size of the ME areas is decreased by a relatively minor 
term (for k B T <~J) which scales with (T - T ) 2 . 

The evaluation of F s (k, u) is detailed in Appendix B. 
In the | a; | e* lin regime, one has F s (k,w) ex w for [4lj] 
T > T« (which is enabled by Eq. This results, 

through Eq. (|B1[) , in the existence of low-energy tails in 
A s (k,Lo), and thus in the disappearance of the spin gap 
and the sharp resonance mode, existing [l6[ at T < Tg, 
in agreement with experiment. 

Typical results for the T dependence of Fg'(k) = 
lim iij ^ [r s (k, w)/w] in the LE and ME areas are pre- 
sented in Fig.ifa). For k B T <C T* a (see Eq. {53}), T s ' 
is characterized by scaling with 1/(T — To) in the LE ar- 
eas, and by a constant plus a minor term cx (T — T)) 2 
in the ME areas (it practically vanishes in the HE ar- 
eas). For k B T ^> TQ a , Tq is characterized by scal- 
ing with 1/[(T - T ) 2 (l - T /T)] in the LE areas, with 
1/[(T-T )(1-T /T)] in the neighboring ME areas (which 
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turn into LE areas at higher T — see Fig. EJa)), and with 
1 /T 2 in higher-energy ME areas. A typical T dependence 
of the k- integrated rg'(k) is also presented in Fig. EJa); 
it is characterized by some linear increase with T at low 
T > Tg, followed by a decrease around k B T ~ T^, and 
a weak dependence at higher T. 

Typical results for r s (k, uj), as a function of uj (at con- 
stant T > Tg), within the LE, ME and HE areas, are 
presented in Fig. EJb) . They are characterized (see Ap- 
pendix B) by scaling with uj for \uj\ <~e s (k) (which is 
somewhat larger than k B T in the LE areas), with a slope 
which is approximately oc l/e s (k), for low e s (k), and 
is smaller for high e s (k); in the case that e s (k) <^J, 
this is followed by a decrease in the \T S \ vs \ui\ slope at 
\lu\ ~ e s (k), and for small e s (k), also by a wide maximum 
of |r s |, centered somewhat above \u>\ ~ e s (k), and its de- 
crease above it; the |T S | vs \ui\ slope rises at e s (k) ~ J, 
and an approximately linear increase of |r s | with \ui\ is 
maintained at higher \uj\ values, up to about the high- 
energy limits of the lagron spectrum (see Fig. [2]). 



5. Non-FL electron scattering rates 



Let e?(k) + iir«(k) and e p (k) + ±ir p (k) be the poles of 
QJ 1 and Q_ s , respectively; they have weight factors w p (k) 
and W p (k) normalized as: J2 p w p (k) = ^ p w p (k) = 1- 
r*(k) and r p (k) are related to the scattering rates 
r 9 (k, uj) and T s (k, w), analyzed above, through the 
Kramers-Kronig relation [39, 40]: 



/ 



dw|r(k,w)|/27r 



(w - ep(k))* + (T p (k)/2)2 



= 1. 



(40) 



This indicates that the dependence of T*(k) on T and 
e p (k) is close to that of r*(k, uj) on T and ui ~ e«(k). 

The elements of the approximate electron Green's- 
function matrix (see Eq. (TK))0 . represented by the 
bubble diagrams in Fig. (Hb), could be expressed as 
MM: 



G d (k,z) = £™ p ,(k> p (k-k')sign(e p ,(k')) 

pp'k' 



A d (k,Uj) 



x cosh 2 (^(e p ,(k') + |z|r p ,(k')|)) 



[b T (e s p ,(k>) + **|r;,(k')|) + / r (-eg(k - kQ - i ? r p (k - k'))] 

[z - (4(k - k) + £p ,(k')) - *i(rg(k - k) + |r p ,(k')|)] 



+ sinh 2 (^( ep ,(k')-^|r p ,(k')|)) 



K(e p ,(k>) - ±*|r p/ (k')|) + fA^(k - k') + *irg(k - k'))] 



(41) 



[z - ( £p (k - k') - e p ,(k')) - Ji(Tg(k - k') + |r p ,(k')|)] 
%g d (k,uj-i0+)/ir = — ^ / dw'A*(k / ,a; / ){A fl (k-k , ,w-w / )cosh 2 (^i,(fa; / )) 

k' ^ 

x [b T (uj 1 ) + / T (a/ - u)] + A q (k - k', u + uj') sinh 2 (& (a/)) [& T (w') + / T (w' + w)]}. (42) 

I 



Since the svivons are in a combination state of 
degenerate condensates (corresponding to the differ- 
ent Q m points), the k' summation in Eq. (j42[) in- 
cludes the averaging of the A s (k',uj') cosh 2 (££, (uj')) and 
A s (k',ui') sinh 2 (££, (uj')) terms over these condensates. 

In the non-FL regime, where major features of the 
physics of the cuprates are described correctly replac- 
ing Q_ d in Eq. (Til?!) by Qi, the evaluation of these fea- 
tures is approached approximating the electron spectral 
functions A d (k,uj) = %g d (kuj - i0+)/ir by A d (k,ui) in 
Eq. (|42|) . Note that the w-integral of A d (k,uj), given in 
Eq. (j42|) . is short of one exactly by the contribution of the 
electron states of the upper Hubbard band, created by 
the second term in the rhs of Eq. Q, ignored in Eq. §5$. 

The evaluation of Ag(k, u>) is detailed in Appendix C, 
and it is found to be, approximately, expressed as a sum 
of two terms, presented in Eq. (IC1|) . One term, Ao c (k, ui), 
represents a convolution of QE and modified svivon spec- 
tral functions, where the low-energy (< k B T) tails of 
A s have been truncated (see Eq. (|B1[) and Fig. Efb)). 



The other term, Ap b (k, uj), represents anomalous effec- 
tive electron bands, generated by the QE spectral func- 
tions and the truncated A s tails. 

As was discussed above, the QE spectral functions 
are specified as low-energy (LE), mid-energy (ME) and 
high-energy (HE) types of functions: Af e (ui), A^ nc (uj) 
and A^ e (ui), respectively; consequently, the A^ h electron 
bands corresponding to them are referred to as LE, ME 
and HE bands. The LE A d h band, due to the arcon 
and stripon peaks, is flat and extends over the low-QE- 
energy BZ areas (see Fig. Efc,f)). The ME A^ h band, 
due to the humpons, merges with the LE band at high T 
(see above) , and has a similar role to that of the HE A^ h 
band at low T. 

In the low \u>\/k B T limit, Aq is dominated (through 
Eq. ([UI])) by the LE A d h band [Uj], and the physical 
properties based on it reflect its anomalous features; thus, 
as is explained in Appendix C, the linewidth of this ef- 
fective band increases linearly with T, and the spectral 
weight W /d (k) within it depends on T (especially in the 



15 



antinodal areas), increasing with it at low T > T', and 
saturating at fc B T > w H ; the band lincwidth plays the 
role of the electron scattering rates r d (k, uj) in physical 
properties derived from it. Typical results for these effec- 
tive T d and W d , as functions of T (for w — >• 0), in k points 
ranging between the nodal and antinodal BZ areas, are 
presented in Figs. [7{a-b). 

On the other hand, in the high-|w|/fc B T limit A d is 
dominated by Aq c (based on a convolution of QE and 
modified svivon spectral functions) [58[ plus a minor con- 
tribution of the ME and HE A d h bands; when T is in- 
creased, the spectral weight within them is decreased by 
the same amount that it is increased within the LE A^ h 
band (see Fig.^b)). Physical properties could then often 
be formulated in terms of QE and svivon contributions, 
where one of them may be dominant. 

Thus, the electron scattering rates T d (k, uj) measured 
in the high-|aj|/fc B T limit, in various physical properties, 
are obtained as convoluted combinations of those of the 
QEs and svivons, as is implied from Eqs. (pIUIFfTj) ; due to 
the cosh 2 and sinh 2 factors in Eq. (HIT) , the major svivon 
contribution comes (see Eq. (|39p ) from LE and ME BZ 
areas at energies ranging between k B T and somewhat 
above e s (k), where a maximum or a plateau appears in 
the V s curves in Fig. Hfb). Consequently (see Figs. (S](b) 
and 13(b)), when T > Tg, T d (k, uj) are approximately 
linear in uj, for k B T < \uj\ < 0.3 eV. 

Typical results for T d (at T above HJ T« and below 

Toa/^B an d wA (Q)Ab)i as a function of uj, in k points 
ranging between the nodal and antinodal BZ areas, are 
presented in Fig. [Tfc). The evaluation of these results is 
based on precise Fermi and Bose distribution functions, 
and not on their asymptotic values assumed in Eq. (jCll) ; 
thus, they reflect a smooth crossover between the high- 
and the low-(|w|//c B T) behaviors, around \uj\ ~ k B T, 
manifested in lower slopes of the T d vs. \uj\ curves for 
|w| ->■(). 

Similarly, due to large cosh 2 and sinh 2 factors, svivons 
of the LE and ME areas also have a major contribution 



to Aq c in Eq. (jCll) . Consequently, the convoluted QE- 
svivon states at energies k B T < \uj\ < 0.3 eV have a band- 
like behavior of linewidth which approximately increases 
linearly with \u\, as can seen in Fig.[3][d-e) [58{ . 

The linear dependence of the linewidths on uj is in 
agreement with ARPES results [13,[5^|; furthermore, the 
measured zero-energy-interpolated linewidths are closest 
to zero along the line of nodes and largest in the antinodal 
direction, in agreement with the results obtained here for 
Tp(k) in nodal and antinodal points (see Eq. (IA3|) and 
the discussion in Appendix A). This linear dependence 
is found here both for the main and the shadow bands, 
as has been observed by ARPES [Uj]. Note, however, 
that Aq misses the effect of the multiple-scattering FL- 



like electron poles in Eq. (|T9]) . which contribute only to 
the main bands (see Fig. E^d)). 

Thus, in the regime where the contribution of the non- 
FL Ag(k, uj) to A d (\s.,ui) provides a good approximation 
for the evaluation of the effective r d (k, uj) in various ex- 
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FIG. 7: Typical results (in arbitrary units), within the non- 
FL regime (corresponding to the T q and V s results presented 
in Figs. [5] and [6}, for the T dependence of (a) the electron 
scattering rates F d , (b) the spectral weights W d within effec- 
tive LE electron band (see discussion in the text), and (c) for 
the uj dependence of F d . 



periments, the measured T d are expected to be charac- 
terized by a linear dependence on T in the low-|w|/fc B T 
limit, and by a linear dependence on uj in the high- 
\ui\/k B T limit. Such a behavior corresponds to the MFL 
phenomenology [Uj], and thus, the hump phase (42| is 
referred to here as the MFL phase. 
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The linear dependence of the r on T or u, down to 
their lowest value within the MFL phase, has been at- 
tributed to quantum criticality [5§, [6(| ■ In the present 
work a different, though not contradictory, approach is 
applied, and the linear dependence on Torw results from 
the role of the Bose svivon field which has linear low- 
energy spectral tails, down to u) = 0, due to Eq. (15B1 . 
This linear dependence, and the switching in the roles of 
k B T and \oj\ between the low- and high-|u;|/£; B T limits, 
result here from the behavior of terms like — |w|6 T (— \w\) 
in these limits. 

The effect of convolution with svivons on the low- 
energy QE spectral peaks, including the different con- 
tributions to Aq, is sketched in Fig. 0|b). The traces of 
the arcon and stripon peaks are considerably broadened 
in the MFL phase [58(; however, since the evaluation of 
physical quantities in the low \uj\/k B T limit is dominated 
by the LE Ag h part of Aq, for which the broadening of the 
QE peaks is ~2k B T (see Appendix C), these quantities 
are sensitive to the sharpness of these peaks. 



6. Anomalous physical properties 

The T dependencies of the transport properties of the 
cuprates, whose anomalous behavior has been associated 
with non-FL behavior from the start, are the consequence 
of the unconventional T behavior of the effective scatter- 
ing rates and spectral weight corresponding to the LE 
Ag b band. The T dependence of the electrical resistivity 
in the cuprates [l[ follows in the MFL phase that of the 
evaluated effective T d , presented in Fig. EJa). 

The effect of transfer of spectral weight from the Aq c 
and ME A^ h bands to the LE A^ h band, when T rises, is 
reflected in Hall-effect and thermoelectric power (TEP) 
results which are dominantly determined by the latter. 
This transfer occurs at both positive and negative ener- 
gies of magnitude ~fc B T; thus, it should be reflected in 
such measurements as an effective increase in the density 
of carriers, subject to the question whether the occupa- 
tion of the LE A^ h band passes through half filling. 

In conventional bands both the Hall number (n H ), and 
the TEP (S), are negative when the band is almost 
empty, and positive when it is almost full; the analytical 
periodic behavior of such bands over the entire BZ gen- 
erally results [3l| in a coincidence between their signs. 
As was discussed above, the LE A^ h band is not con- 
ventional, and consists of the arcon-stripon contribution 
A? (lj) to the QE spectral functions (see Figs. 1314] ). and 
the low-energy tails of the svivon spectral functions (see 
Fig. [6]). Thus, it occupies only a part of the BZ, and a 
careful analysis should be made on the origin of the signs 
of n H and S, derived from it. 

The sign of S is [6l[ that of the average energy (rel- 
ative to the chemical potential) within the range of the 
band close to the FS, introduced through ujdf T (oj)/dui. 
Since the LE A^ h band lies almost entirely within this 
range, the TEP derived from it could be either positive 



or negative, depending on the band occupation, and thus 
on stoichiometry. 

On the other hand, the sign of n H is [6l[ opposite 
from that of the averaged second derivatives of the band 
dispersion curve over the FS (where —df T (uj)/duj con- 
tributes). Since the LE A^ h band extends only over the 
arcon-stripon QE BZ areas, shown in Fig.[3ja-c), where 
the sign of the second derivatives of the band dispersion 
curve is dominantly negative, the sign of n H is expected 
to be positive, independently of the occupation of the 
band (and thus of stoichiometry); the fact that this low- 
energy band is flat results in small dispersion derivatives, 
but also in a nonzero —df T (u>)/duj factor over a range of 
the BZ, and the contradicting effects introduced to the 
expression 61] for n H compensate for each other. And 
indeed, for hole-doped cuprates, studied here, n H is found 
to be positive in the MFL phase, and its T dependence 
follows @,[3 that of W d {k), presented in Fig.^b). 

The T dependence of the TEP in YBCO has been ana- 
lyzed Q in terms of a "narrow-band model" , under which 
IS* | is larger than typical metallic values, and increases 
with T up to its saturation value (when k B T exceeds the 
bandwidth), thus [4|: Ssat = {k B /e) In [x/(l — x)]; here x 
is the fractional band occupation (thus x — when the 
band is empty of electrons, and x = 1 when it is full) 
within the measured stoichiometry. 

For hole-doped cuprates of an electronic structure sim- 
ilar to the one presented in Fig. [31 the TEP does not sat- 
urate with T , and has an almost universal dependence on 
T and the doping level [5|. Thus, in the underdoped and 
lightly overdoped regimes, S increases with T at low T, 
reaching a positive maximum, and then decreases with T, 
having an almost linear dependence on it at high T; at 
the optimal stoichiometry, S crosses zero at T = 300 K 
which has been applied [5J to determine the doping level 
in cuprates. In the heavily overdoped regime, S is neg- 
ative and decreases monotonously with T, consistently 
with an FL state with a band which is less than half 
filled (as is expected in this regime). 

The LE A$ h band is a narrow one (being based on 
the almost dispersionless stripons and arcons); however, 
since its linewidth increases linearly with T, the TEP sat- 
uration temperature [4[ cannot be reached. Furthermore, 
the asymmetry of the QE and electronic spectra, pre- 
sented in Fig. [3J with respect to the inversion of the sign 
of the energy, results in a temperature-induced transfer 
of spectral weight to the LE A^ h band (from the Aq c and 
ME A^ h bands) which is larger for negative than for pos- 
itive energies. Consequently, the fractional occupation x 
of this band decreases when T is increased. 

Thus, in the underdoped and lightly overdoped 
regimes, the LE Ai h band is more than half filled at low 
T , resulting in positive S (within the MFL phase) which 
first increases with T to approach its saturation value, 
but the decrease in x with T results also in a decrease in 
this saturation value Q , and thus in S (after it reaches its 
maximal value). Consequently, a "universal" behavior of 
S is obtained for cuprates of the same type of electronic 
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structure as in Fig. |3l in agreement with experiment Q. 
The value S = 0, at T = 300 K, at optimal stoichiometry, 
implies that the LE A°] h band is then half filled. 

The existence of a universal behavior in the cuprates 
has been observed @ in various physical quantities, de- 
termined by the LE A d h band. This behavior is expressed 
through the scaling of their T-dependencies in a typi- 
cal energy corresponding to spectrum of the low-energy 
QEs. These quantities include, in addition to the TEP, 
the planar and c-axis resistivities, the Hall coefficient, 
the uniform magnetic susceptibility, and the spin-lattice 
relaxation rate. 

Linearity of the electron scattering rates in uj, for 
|cj| < 0.3 eV (similarly to ones shown in Fig. Etc)), 
within the MFL phase, has been observed in optical re- 
sults [iiilsH- In different hole-doped cuprates, the op- 
tical conductivity a d (uj) consists in this phase [T3| of a 
low-energy Drude term, and a higher-energy "mid-IR" 
term, where the effective density of carriers within the 
Drude term is 4-5 times smaller than those integrated 
up to above the mid-IR energy (a behavior referred to as 
"Tanner's law"). 

Within the present approach, the electron spectral 
functions are approximated by Aq(\c,ui) in Eq. ((^.in- 
cluding both the Ay h and A d c terms in Eq. (IC1|) [58j . 
At low w, the QE contribution to A d includes only the 
stripons and arcons peaks, resulting in the Drude term 
in a d {uj). The mid-IR term in a d {uj) is introduced due 
to the growing role (in Aq) of humpons at higher ener- 
gies, and of the almost detached higher-energy QEs (see 
Fig.[2Ia-b)) at further higher energies. Thus the increase 
in the effective number of carriers with energy, expressed 
by Tanner's law [l(|, reflects the difference between the 
occupied QE spectral weight residing within their low- 
energy peaks, and that residing within their entire spec- 
trum. 



D. Paired QEs and electrons 

1. Formulation and general features 

QE pairing is approached adapting of the Migdal- 
Eliashberg theory ,27) for the present case. Nambu 
spinors (?^(k), q (— k)) are considered as creation oper- 
ators (the k values here are within a half of the BZ, but 
since there is a freedom in choosing this half, the derived 
expressions are valid within the entire BZ). Self-energy 
2x2 matrices are derived, and expressed (in terms of the 
Pauli matrices r 1 , r 2 and r 3 ) through: 

eg(k) + E«(k,*) = [ E g(k) + E'(k,z)]r 3 + »$»(k,«)T 1 
+3**(k,z)r 2 . (43) 

The diagonalization of ejj(k) + E 9 (k, z) yields QE Bo- 
goliubov states annihilated by <7 + (k, z) and q_(k,z), 
corresponding to positive and negative energies, re- 
spectively; the transformation between them and the 



unpaired-state QE operators is expressed as: 

q(k) = exp^(z)/2)[cos(^(z))q + (k,z) 
-sin (£(z))q_(k, z)], 
gt(-k) = cxpH$(z)/2)[sin(£(z))g + (k,z) 

+ cm (k,*)]. (44) 

This transformation diagonalizes Dyson's equation [39|, 
EoT ]. yielding the poles of the diagonalized QE Green's 
function Q±(k,z) at: 

z = ±e«(k) + E^(k,z) 

= ±E q (k,z) + iSE^(k,z), where 
QE^(k,z) = ±cos(2££(.s))3E 9 (k,z), 



E q (k,z) = v /^(k)+^(k,z)]2 + [$9(k,z)]2, 
$«(k,z) = |$«(k,2)|exp(*V£(aO), (45) 

K( z ) = ^l( z ) or i>w( z ) + 7r ^ and 
$«(k,z) = cos(4>l(z)-cf>l(z))\<P q (k,z)\ =±|$ 9 (k,z)|. 

The coefficients ££(z) are obtained through: 



cos(2^(z)) 



e q {k) + KZ q (k,z) 



sin(2^(z)) 



$«(k,z) 
S9(k, z)' 



(46) 



The spectral functions of the QE Bogoliubov states are 
obtained through: 



A q ± (k,w) ee %g q ± (k,w-i0+)/ir 
± 



r^(k,w)/27r 



[ W T^(k,^)] 2 + [^r 9 ± (k, w )] : 



(47) 



where T q ±(k,uj) = 23E^(k,w-i0+) > 0, and thus A q ± > 
(note that cos(2^(o;)) is positive within the uj > 
range of A q + , and negative within the uj < range of 
A q _— see Eqs. gSEU)). 

The self-energy terms E 9 and $ 9 are primarily deter- 
mined, at low energies, by QE-lagron coupling, through 
AH in Eq. (TTS|) ■ Their self-consistent expressions are de- 
rived similarly to the un-paired case in Eq. (|34p . yielding: 

£"(k,z) - [ d^[cos 2 (^_ q K)) 

q J 

xA 9 + (k-q,^)+sin 2 (^_ q (^)) 
xAi(k-q,^)]^ r (q,^,z), 

<F(k,z) SS ±-J2 f dw«frin(2#_ q ( W ')) 
q J 

xexp(i^_ q ( W «))[^.(k-q, W «) 
-A q _(k- q ,Lo q )]F T (ci,oj q ,z). (48) 

5E 9 (k, uj — i0 + ) is obtained by expressing (uj — i0 + =F 
w A (q) -uji)- 1 in F T (q,w«,w-i0+), in Eqs. ([M1I3S]). as 
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[39l lioj 7riJ((cj =p w A (q) — w 9 ); such an expression could 
not be applied to obtain $5<I> 9 (k, ui — i0 + ), because of the 
exp (i0k_ q (w 9 )) phase factors in the rhs of Eq. (|48]) . but 

it could still be applied to obtain a term i$ 9 (k, uj — i0 + ) in 
$ 9 (k,w-iO+). Due to the f T and 6 T factors in Eq. (p4|). 
QEs and svivons at energies w 9 and w A (q) contribute to 
3£ 9 (k, u - i0+) and $ 9 (k, u - i0+) at w = uj q ± w A (q) 
whenw A (q) <C k B T, and at w = w 9 +sign(u; 9 )u; A (q) when 
a/(q)»fc B T. 

The most significant low-energy contribution to £ 9 and 
$ 9 in Eq. ([15]) is from q points corresponding to Q m and 
Q m -vicinity lagrons (see above) ; an additional significant 
contribution comes from Q-ESP lagrons. Thus, by Fig. [2] 
a QE spectral peak around the energy uj q , located at ~k— 
Q m , contributes (due to Q m and Q m - vicinity lagrons) 
peaks to 3E 9 (k, u - i0+) and $ 9 (k, u - i0+) atw- w q , 
with asymmetry towards larger |w 9 |; if such a QE peak 
is located at ~k Q, it contributes to them (due to 
Q-ESP lagrons) peaks at uj ~ uj q + sign(w 9 )u; A (Q), if 
k B T < w A (Q), and at u ~ w 9 ±w A (Q), if fc B T » w A (Q). 

SRE 9 and <1> 9 could be expressed through the Kramers- 
Kronig relations [1^, H(| : 



5ft£ 9 (k,w) 



$ 9 (k, W -zO+) = p / =1^. 



dw'3E 9 (k,w' -i0+) 

it(lu — tjj') 

dw'$ 9 (k,w' -i0+) 



(w - w') 
+i$ 9 (k,w-i0+). 



(49) 



The effect of peaks in 3X 9 (k, w — i0 + ) (which is always 
positive) on 5ft£ 9 , in Eq. (j49|) . is to "push" the peaks in 
-A±(k, lu) away from them, through Eqs. (|4"5l |4"T1) . As was 

discussed above, the location of peaks in $ 9 (k, ui — i0 + ) 
(and in 3I] 9 (k, u — i0 + )) is close to k — Q, and sym- 
metrically around it. Consequently, in order to stabilize 
pairing, one should choose (see Eqs. (1451) . (l46l) and (|4"8]l): 

0k M = ^k( w )j getting 
^k+Q ( w ) = 0k ( w ) + 7r ! an( i introduce 
<l 9 (k,w) = |$«(k,w)|cos<^(u;), where 
sign[<l 9 (k, a;)] is independent of uj. (50) 

$ 9 (k, w) plays the role of the pairing order parameter; 
by Eq. (|50"|) it reverses its sign when k is shifted by Q, 
resulting to an approximate d x 2_ y i pairing symmetry. 
For the choices in Eq. (1501) one gets that the effect of the 
peaks in <l 9 (k, u - i0 + ) on $ 9 in Eq. ((49]) is to "push" 
(through Eqs. (|45[|4"T)l ) the peaks in A!j_(k, uj) away from 
them, if they are at energies of the opposite sign, and to 
"pull" these peaks towards them, if they are at energies 
of the same sign. 

Such a scenario could result in a pairing gap, where in 
addition to the humpons, the stripon and arcon peaks 
(which are not on the lines of nodes) are split into 
positive- and negative-energy peaks. As was discussed 
above, the locations (at low T) of the centers of the peaks 



in 3£ 9 (k, uj - i0 + ) and $ 9 (k,u; - i0 + ) are at higher val- 
ues of \uj\ than those of such split peaks in ^4±(k — q> w )> 
where q = Q or Q m . 

By the symmetry of the QE spectrum under a shift in 
Q, a split peak in A q + (\t,uj) or A q _(k,uj) is pushed, on 
the average, away from zero by the opposite-energy-sign 
peaks in 5I] 9 (k, uj — i0 + ), and towards zero by the same- 
energy-sign peaks in it, and this is insufficient to stabilize 
the gap, unless the symmetry between the degenerate la- 
gron and svivon condensates (see Fig. [21 Eq. ([201) . and the 
above discussion) is broken, and a static striped structure 
sets in. 

The additional effect of peaks in $ 9 (k, uj — i0 + ), at both 
positive and negative energies, is to push this split peak, 
on the average, away from zero. Thus, a gap is stabilized, 
under a certain temperature (in a state of combined la- 
gron and svivon condensates corresponding to fluctuating 
inhomogeneities), by breaking symmetry through the in- 
troduction of a nonzero $ 9 . A quantitative calculation 
on this pairing scenario will be published in a separate 
paper. 

Nonzero <3E> 9 results in a nonzero anomalous QE Green's 
function matrix T q (see above), and consequently in 
a nonzero anomalous electron Green's function matrix 
T_ . The expression for its zeroth-order term Tq, in 
terms of T s and , is presented diagrammatically in 
Fig. [TJb). Within the k representation, the expression 
for .7o(k, z) is similar to the expression for Q^{k. : z) in 
Eq. (|4"Tj) ; the cosh 2 and sinh 2 (££, ) factors there are 
replaced by sinh(2££,)/2, and additional factors are in- 
cluded (for both F$ and Go) due to QE pairing; thus, 
Qq (k, z) and ^(k, z) are expressed in terms of the di- 
agonalized 0|.(k, z), and have their poles specified in 
Eq. (05]). 

As will be detailed elsewhere, the electron spectrum is 
determined through the diagonalization of 2 x 2 matri- 
ces whose diagonal terms are determined by Q d (k, z)" 1 , 
and non-diagonal terms by J- d (k, z)~ x and [J" d (k, z)^ 1 }* . 
The evaluation of the approximate electron spectral func- 
tions Ag(k, to) in a paired state is analogous to their eval- 
uation, through Eqs. (14*21 IC1[) . in the unpaired case; but 
the cosh 2 an d sinh 2 (££') factors there are modified, 
and expansion coefficients due to the pairing diagonal- 
ization procedures are included. The obtained electron 
spectrum is characterized by a pairing gap of the same 
symmetry as that of <3E> 9 ; the opposite signs of <J> 9 around 
the lines of modes result in a zero electron pairing gap on 
these lines, and in its continuous variation around them, 
due to QE-svivon convolution. 



2. Pseudogap phase 

As was discussed above, the number of poles in £/ 9 (k, z) 
is different for QEs in nodal and antinodal k points. Thus 
these two types of low-energy QEs correspond to differ- 
ent symmetries, and it is analytically possible to have 
<t q 7^ for only one of them (corresponding to a specific 
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area within the QE BZ). Since (as was discussed above 
with respect to humpon-humpon splitting) Q m lagrons 
contribute significantly (through Eq. ([4*8)1 ) to QE pairing 
in antinodal points, and less (if at all) in nodal points, it is 
likely that there exists a temperature range T c < T < T* , 
where $ 9 ^ in the antinodal, but not in the nodal QE 
BZ areas. The symmetry of a state within this T range is 
different from the symmetries of the states where <I> 9 ^ 
for either almost all the low-energy QEs, or almost none 
of them. 

Scattering between the <I> 9 ^ antinodal, and the 
$ 9 = nodal QEs prevents the existence of supercur- 
rent in such a partial-pairing state, and it corresponds 
to the PG phase of the cuprates [42|. From a compari- 
son between the expressions for £ 9 and $ 9 in Eq. (jUJ), 
and the above discussion, one concludes that T* could 
not exceed a maximal value which is somewhat smaller 
than ^uj H /k B (where the effect of the inhomogcncitics 
disappears). Consequently, k B T* could be, at the most, 
close to ~J (see discussion above), in agreement with the 
values of T* observed in the phase diagram [42| . 

The stripon peaks split in this phase into positive- and 
negative-energy peaks, as is sketched in Fig. 0Ja). The 
contribution of Q-ESP lagrons to and $ q in Eq. gSJ), 
and the fact that |5ft£ 9 | > |$ 9 | within the higher- weight 
Bogoliubov band, results in an energy separation [56] > 
w A (Q) ~ jE lea between same-energy-sign stripon and 
humpon peaks. 

There exists a low-T regime within the PG phase where 
the inhomogeneities become static (resulting in an en- 
ergy gain), and they are observed as a glassy structure 
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]; this regime is characterized by a lower mini- 
mum e^ in of the svivon dispersion curves, and a higher 
spectral weight in the humpon peaks (which, as was dis- 
cussed above, correspond to the inhomogeneities), on the 
expense of that in the split stripon peaks. In order for a 
static inhomogeneous structure to coexist with pairing, 
it should not break the symmetry between the four Q m 
points in Eq. (|2"0j) and Fig. [*"] and e^ in should not be too 
close to zero; this requires the effect of coupling to the 
lattice, manifested in the existence of local distortions. 

The static inhomogeneities are obtained through two 
types of combinations of striped condensates, corre- 
sponding to different Q m points; one type is of the differ- 
ent phases of the SDWs introduced by the stripes, result- 
ing in static charge-density waves [22| (CDWs) of wave 
vectors 2Sq m (while the spins are fluctuating); the other 
type is of the different directions (thus a and b) of the 
stripes, resulting in a checkerboard-like structure [2l[ of 
short stripe segments of ~4a periodicity, along the two 
directions (similarly to a structure predicted earlier 



18 



35]) 



Longer unidirectional stripes can be obtained 
when the resulting energy gain is larger than that due 
to pairing which is suppressed altogether, as probably 
occurs in "1/8 doping anomalies" and below the minimal 
SC doping level. 

The above derivation of the 1ow-cj QE and svivon scat- 
tering rates in the MFL phase required the existence of 



low-w QE (particularly stripon) spectral weight, as oc- 
curs [41| for T > Tg. Consequently, the opening of a 
partial QE gap in the PG phase results in a substantial 
reduction in T q and |r s |, for below the gap energy, 
but not above it, in agreement with experiment [HQ. 

However, since the gap is partial, the low-w scattering 
features, including the low-energy svivon spectral tails, 
existing for T > T* , partially persist for T c < T < T*. 
Consequently, the electron spectral functions in the PG 
phase include, as in the MFL phase (see Eq. (|C1|I ). an 
effective LE band part ^^(k, w) which consists here of 
the low-energy part of the Fermi arcs around the lines of 
nodes, and a convoluted QE-svivon part AQ C (k,uj). The 
contribution of arcons to Aq c forms the extension of the 
Fermi arcs to higher energies, and the contribution of 
stripons to it is modified by pairing (see above), reflect- 
ing the existence of the pseudogap. The association of 
the latter electron states with Bogoliubov states, due to 
pairing, is supported by experiment [63l. l64j . Also the ex- 
istence of Fermi arcs which are distinct from Fermi pock- 
ets (see above) has been established by experiment [65$ . 
These two features in A d in the PG phase are sketched 
in Fig.^b). 

A decrease in T, within the PG phase, results in the 
transition of the spectral functions of low-energy QEs 
from a structure of $ 9 = arcons to that of $ 9 ^ 
stripons and humpons. It is associated with the grow- 
ing effect of the inhomogeneities at lower T, within the 
PG phase, discussed above. Since all the arcons are cou- 
pled to stripon-humpons through lagrons, the transition 
of unpaired arcons into paired stripon-humpons would 
include at T — all the arcons (except for those on the 
lines of nodes, where $ 9 = by symmetry) if the arcons 
remained unpaired. 

The signature of such a transition on the electron spec- 
trum is the observed reduction of the Fermi arcs with 
decreasing T, within the PG phase, and indeed, the ex- 
trapolation of their reduction to T = results in the 
points on the lines of nodes (66[. Such a low-T "d-wave 
nodal liquid" has been observed [67j in the regime where 
the PG state persists down to T — > 0. 



3. Superconducting phase 

SC occurs when <J> 9 ^ for the nodal arcons (except 
for those on the lines of nodes) , in addition to the antin- 
odal stripons and humpons. Thus the arcon peaks split 
below T c into positive- and negative-energy peaks, as is 
sketched in Fig. 0Ja). Since their coupling to other QEs 
through Q m lagrons is weak, or absent, the additional 
pairing necessary for SC to set in is induced primarily 
through Q-ESP lagrons. 

As was discussed above, the lagron energy involved 
in such pairing is ^w A (Q) which could be treated analo- 
gously to the dominant boson energy [68] k B fi in electron- 
phonon systems. A coupling parameter A is obtained by 
approximating Eq. (|48p similarly to the treatment of such 
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systems [68|, resulting in values as large as A ~ 3. The 
dependence of A on the doping level is weak since, below 
T e , $ 9 7^ both for the nodal and the antinodal QEs. 

Estimating T c on the basis of an expression (for 
electron-phonon systems) which is valid within the en- 
tire coupling regime [6§|: T„ ^ 0.25fi/[exp (2/A) - l] 1 / 2 , 
results in an approximate scaling w A (Q)/Ar B T c — 3.9, and 
thus [56[ E iBB /k B T c ~ 5.2, in agreement with experiment 
[69| . This explains the enhancement of T c compared to 
conventional electron-phonon SCs, where such large A 
values would result in lattice instabilities. 

The crucial role played by lagrons at the Q m points 
and their vicinity (where w A (q) oc |q— Q m |) in the pairing 
of stripons is not comparable to the minor role played by 
low-energy acoustic phonons in conventional SCs; while 
the coupling constants between electrons and acoustic 
phonons vanish [3lll39| for q — > 0, the coupling constants 
7(q) between QEs and lagrons remain constant for q — > 
Q m (see Eq. (IT5)l and the above discussion). 

High-energy lagrons also play some role in pairing, 
and specifically the high-energy SP lagrons (see Fig. 
An optical analysis, in an attempt to determine the en- 
ergy spectrum of the bosons involved in pairing in the 
cuprates, reveals [70l ] significant contributions of bosons 
at low energies, at energies [56J ^w A (Q) ~ jE ree , and at 
energies ~0.2-0.3 eV which are consistent with those of 
high-energy SP lagrons. 

The opening of an SC gap prevents the above deriva- 
tion of QE and svivon scattering rates, depending on 
the condition T > (which is only partially fulfilled 
in the PG phase). This results in a drastic reduction in 
T q and |r s |, for \uj\ below the gap energy, but not above 
it. Consequently, the split stripon and arcon peaks are 
sharp below T c , while the humpon peaks remain wide 
(see Fig. mja)). 

Also the low-energy svivon spectral tails disappear be- 
low T c , and thus the electron spectral functions A d (k, u), 
can then be approximated through a modification (see 
above) of the convoluted QE-svivon term Ag c (k, w) in 
Eq. (jCll) alone. The resulting radical reduction (for 
T < T c ) in the electron scattering rates, for \lu\ below 
the gap energy, is observed in optical @ and microwave 
results. 

An aspect of Tanner's law discussed above, is 

that the optically determined SC-state superfluid den- 
sity corresponds to the density of carriers determined by 
the normal-state Drude term, and misses the contribu- 
tion of the mid-IR term. This observation is consistent 
with the present prediction that both the Drude and the 
superfluid density are determined by the contribution of 
stripons and arcons to A$. 

The traces of the arcon and stripon gaps on A d , 
sketched in Fig. H£b), are often referred to as nodal 
and antinodal gaps, respectively. The antinodal gap- 
edge structure includes the (stripon-derived) peak and 
(humpon-derived) hump (see Fig. BJb)). The nodal gap 
has been found [zHlzH to scale with ~5fc B T c , consistently 
[68j with the ratio between the pairing gap and T c for the 



large coupling parameters A in this case (see above). 

Since Eq. P5|) couples between the nodal and antin- 
odal pairing order parameters, the opening of a nodal 
gap below T c modifies the antinodal gap (which opens 
below T*) in a manner observed (73l. l74j as an apparent 
onset (below T c ) of a "second energy gap" , superimposed 
on the pseudogap; however, it is clear that the antinodal 
gap is "smoothly connected" to the nodal gap [75J. 

In the underdoped regime, most of the antinodal (but 
not nodal) Cooper pairs are formed at T c < T < T*; 
this is consistent with the observation [76| that in this 
regime only a small density of additional antinodal (but 
not nodal) pairs are formed at T < T c . Since the super- 
fluid density consists of pairs formed both at T < T c and 
at T c < T < T* , it does not coincide with the density of 
the additional pairs formed at T < T c [76j . 

The svivon energy minima are lower below T c than 
their values (e^in) above it, and since the low-energy 
svivon spectral tails existing above T c (see Fig.|6]and Ap- 
pendix B) have disappeared, a spin gap opens up. This 
results 16] in the existence of a sharp resonance mode be- 
low T c . Since (see above) the energy separation between 



the stripon and humpon peaks 



> ^(Q) ~ *E„ 



the reduction in the width of the stripon peak results in 
the appearance of a dip between the peak and the hump 
(see Fig. ffl^b)), at ~E above the peak, in agreement 
with experiment [H, [z3] ■ This structure is manifested as 
a T-dependent "antinodal kink" in the electron disper- 
sion curves 



78] 



The balance between energy gains due to pairing and 
inhomogeneity, discussed above with respect to the PG 
phase, is relevant also for the SC state, especially in the 
underdoped regime. The fact that the lattice is involved 
in the establishment of the static checkerboard-like [2l[ 
structure (see above) is supported by the observation of 
an unconventional isotope effect within this regime [79| . 
It is viewed not only in T c , but also in properties like 
the magnetic penetration depth, in T*, and in the onset 
temperatures of spin-glass and AF ordering in adjacent 
inhomogeneous regimes of the phase diagram p79j . 

Since the local lattice distortions have a non-linear 
negative effect on the pairing energy, a heterogeneous 
nanoscale structure sets in [80( , including regions of vary- 
ing strengths of the effects of pairing and inhomogeneity. 
Regions of a stronger pairing effect are characterized by a 
larger spectral weight within the peak (on the expense of 
that within the hump), and a weaker checkerboard- like 
structure. Since the static inhomogeneities have only a 
minor effect on the nodal gap, this heterogeneous struc- 
ture has almost no effect on it (it has a significant effect 
on the antinodal gap), in agreement with experiment [80| . 



V. CONSTRAINT SUSCEPTIBILITY 

The spectrum w A (q) of the lagrons, and the constants 
7(q) of their coupling to QEs and svivons, have to be 
determined through the condition that the resulting QE 
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and svivon spectra satisfy, in every site, the auxiliary- 
particles' constraint in Eq. ©. This condition can be 
expressed as a requirement [la] that two-site correlation 
functions of the svivon operators in the lhs of Eq. ([5]) is 
equal to that defined, similarlythrough the QE operators 
in the rhs of Eq. ([6]); thus: 

^2( s L s it7 a j<r> s ja') ~ (QiQi9jQj)- ( 51 ) 

go' 

The terms in the lhs and the rhs of Eq. (|5Tj) introduce a 
susceptibility-like function, referred to as the constraint- 
susceptibility [l6[ (Xc or Xc); h should be the same (thus 
Xc = X?)> whether it is evaluated on the basis of the 
svivon spectrum (thus Xc), or the QE spectrum (thus 
Xc), through the lhs or the rhs of Eq. (|5Tj) . respectively. 
Vertex corrections to Xc an d Xc, due to AH in Eq. (TT2")) 
and W in Eq. (|23p . do not vary them qualitatively. 

A major feature of QXcCq, u ~ in the SC state, is 
fl6j the existence, around q = 0, of a sharp low-energy 
peak at w ~ E zec , and a higher-energy tail, extending up 
to ^J; the Xc = Xc equality implies that such a peak 
structure must be also a major feature of ^xKl, u ~ 
i0 + ) in the SC state; it can be, approximately, evaluated 
through [HUSH (see Eqs. (|33]). <|35]l. (|I7]> and USTJl): 

X q M,z) = / / ^ 2 {[cos 2 (^K))^ + (k, Wl ) 

+ sin 2 (^^ 9 ))^_(k, Wl )][cos 2 (^_ q (c 2 )) 

x^(k-q, W2 )+sin 2 (^_ q ( W2 )) 

xA 9 _(k-q, W2 )] + isin(2^( Wl )) 

x sin (2^_ q («a)) cos [#( Wl ) - #_ q (wa)] 

x L4 9 + (k, Wl ) - AI (k, wi)] [A« (k - q, wa) 

-A«(k-q,u, 2 )])M^AM. (52) 
J z + w 2 — Wl 

The structure of a peak in S?x<;(q = 0,w — «0 + ), ob- 
tained through Eq. (f52"j). represents an average on the 
same-k-point transitions between negative- and positive- 
energy QE spectral features (thus on the two sides of 
the SC gap, as is sketched in Fig. BJa)). These features 
consist of single sharp arcon peaks in the nodal areas, 
and of separate stripon and (broad) humpon peaks in 
the antinodal areas. 

Thus, an equality between 9x1 an d ^Xc, at q = 
in the SC state, means that the sharp low-energy peak 
corresponds to averaged same-k transitions between the 
positive- and negative-energy arcon peaks, while the 
higher-energy tail corresponds to averaged same-k transi- 
tions between the positive- and negative-energy stripon- 
humpon spectral features. The proximity of the low- 
energy peak to E reB , and of the extent of the higher- 
energy tail to ~ J, is consistent with the above-mentioned 
values of the SC nodal gap [11], and the humpon en- 
ergies, respectively. This demonstrates the qualitative 
self-consistency of the theoretical approach, the lagron 
spectrum, and the spectral functions presented here. 



The constraint-susceptibility analysis demonstrates 
that the connection between features (such as typical en- 
ergies) of the QE and lagron spectra, derived here, is 
closely related to the auxiliary-particles' constraint. This 
results in similarities between spectroscopic features de- 
rived here and in models where electrons are coupled to 
spin fluctuations. However, since such models do not cor- 
respond to the large-Z7 case (studied here), a quantita- 
tive analysis based on them [81| does not yield a realistic 
spectrum of the relevant spin fluctuations. 



VI. CONCLUSIONS 

In conclusion, a theory for highly correlated layered 
SCs, where the low-energy excitations are approached in 
terms of combinations of atomic-like electron configura- 
tions, rather than approximately independent electrons, 
has been shown to resolve qualitative mysteries of the 
cuprates. A Lagrange Bose field which accounts for the 
tendency of the doped system to form stripe- like inhomo- 
geneities, enables treating these configurations as bosons 
or fermions. 

The addressed anomalous properties of the hole-doped 
cuprates include the observed phase diagram, non-FL to 
FL crossover, the existence of MFL critical behavior and 
a PC phase with Fermi arcs, kink- and waterfall-like spec- 
tral features, the drop in the scattering rates in the PG 
phase, and further in the SC phase, an effective increase 
in the density of carriers with T and w, the correspon- 
dence between T c , E rcB , and the SC nodal gap, etc. 

Electron-lattice coupling is not included in the Hamil- 
tonian treated here; however, such coupling is neces- 
sary to explain the establishment of static inhomogeneous 
structures within the PG and SC states. This coupling 
is expected to be strong for the derived low-energy spec- 
tral peaks, introducing phonon anomalies at comparable 
energies, an anomalous isotope effect [79], etc. 

The electronic structure of the FeSCs differs from that 
of the cuprates. However, the similarity between their 
anomalous properties [l6j |. and the fact that a formally 
common many-body Hamiltonian could be worked out 
16] for both systems, implies that part of the conclu- 
sions drawn above about the physics of the cuprates, ap- 
ply also for the FeSCs, with some modification. Further 
details of the theory, and a comparison between its con- 
sequences for the cuprates and the FeSCs, will appear in 
forthcoming papers. 
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Appendix A: Evaluation of the QE Scattering Rates 

The behavior of T q is determined at k points within the 
low-energy nodal and antinodal areas, and the adjacent 
high-energy areas; it could be expressed as a sum of two 
terms: 



r«(k )W ) = r;(k, 



r|(k,o;), 



(Al) 



corresponding to the contributions of the Q m lagrons 
(r^), and of the other lagrons (T q ), to the q summation 
inEq. {35]). 

By Eqs. (|2D1I55]>. r?(k,w) in Eq. (JXTJ can be approx- 
imately expressed as: 



r«(k,w)<x£;A9(k-Q m ,w), 



(A2) 



m— 1 



where the combination in the rhs is of terms specified as 
A 9 e (ui), A^ e (w) and A^uj) (see above), consisting of the 
arcon, stripon and humpon peaks, shown in Fig. Eta), 
and the higher-energy spectrum, shown in Fig. G2a-b); 
terms scaling with them contribute, through Eq. (|A2I) . 
to low-|o;| T q within the antinodal areas, and partially 
within the nodal areas (not too close to the lines of nodes) 
and the adjacent high-energy areas. A 9 c (uj — > 0) (and its 
contribution to r\(uj — > 0)) decreases when T increases 
(see below), and lim T ^ A q c {uj — > 0) does not vanish (for 
unpaired QEs). Consequently, in k points within the 
above areas there exists a nonzero limit: 



rg(k) = lim r?(k,w -¥ 0). 



By Eq. (|32|) . and the lack of contributions due to T% 
(see below), T^k) determine the zero-T limit of the 
widths of the stripon and arcon peaks which are, con- 
sequently, ~0 close to the lines of nodes. Let be a 
typical value of T^k) in the antinodal areas; since it is de- 
termined by the coupling of QEs, through Eqs. (|35IA3I) . 
to Q m lagrons, while uj h is determined by their coupling, 
through Eq. ([M]) . to both Q TO and other lagrons (see be- 
low) , Tg a is few times smaller than uj h . 

A k-dependent low-energy scale w L (k) exists which is 
close to Tg a in the antinodal areas; for k B T & \w\ <C 
w L (k), A q (k, uj) is dominated by the contribution of 
Fp(k) to Eq. (|32|) . and it is almost independent of T 
and uj within this range. For k B T 3> w L (k), the T- and 
w-dependencies of the stripon/arcon peak center (dis- 
cussed above) , and of T q , largely determine the behavior 
of A q e (uj), through Eq. ff32]) . Its low-w T-dependence is 
determined by the broadening of the peak (through T q ) 
which, self-consistently, scales with T (see below), result- 
ing in an approximate scaling of the low-o; A q c with 1/T. 

r|(k, uj) in Eq. (|A1[) includes the part in the rhs of 
Eq. (P_5j where only the four Q m points in Eq. (|20p are 



omitted from the sum over q. This summation requires a 
2D integration which could be carried out by dividing the 
BZ into a mesh based on qii and q_i_ lines, directed along 
lagron-energy gradients, and perpendicular to them, re- 
spectively. 

As can be viewed in Figs. [2] and Etc), for a sufficiently 
fine mesh, each BZ section (specified as s) within it could 
be approximated as a circular segment where uj x (q) is lin- 
ear in q|i , and A q (k — q, uj =f w A (q)) in Eq. (|35|) could be 
expressed as A q (uj^fuj x (q)), where A q is of the form of ei- 
ther Af e or A q no or A^ Q . Consequently, the q integral over 
the circular segment s is, approximately, proportional to 
an integral on duj x , where the integrand is multiplied by 



■ Co', 



), between positive limits uj Xq and uj Xq (wp 9 



Us 



a constant) 



Fig. [^indicates that 



M 

J 0s 



UJ 



„? = within four circular 



BZ sections around the V-shape lagron energy minima 
at the Q m points, and that \uj^ q \ are large in integration 
sections corresponding to Q-ESP lagrons (where uj x ~ 
u; A (Q)) and to high-energy SP lagrons. 

The major features of T^k, w), in the low- and high- 
\uj\/k B T limits, are obtained when the following approx- 
imations are applied within each of the above duj x inte- 
gration ranges: (a) the asymptotic behavior of b T (oj x ), 
for uj x < k B T and uj x > k B T, (thus b T {uj x ) = k B T/uj x , 
for uj x < k B T, and b T (uj x ) = 0, for uj x > k B T) is ex- 
tended to uj x < k B T and uj x > k B T, respectively; (b) the 
low-T limit of f T (u') (thus f T (co') = 1, for uj' < 0, and 
f T {uj') = 0, for uj' > 0) is applied for f T (uj x T w). This 
yields, through Eq. (f35|) : 



r«(k,w) = ^[r^ s (k,a;)+r« 8 (k )W )], where 



(A3) r^ s (k, W ) cx k B T / " S duj x (l-uj* q /uj x )[A q s (uj-uj x ) 



+A q Juj + uj x )}, 



(A4) 



r* B (k, w ) 



duj x (uj x -uj* q )A q (uj-uj x sign(uj)), 



M _ 

J Tns - 



min (k B T, uj Xq ), uj^ls = min (fc B T, uj Xq ), 



\ q _ 



= niin(| W |, W A s 9 ), uj Xq xs = min(\uj\,uj Xq ). 

By the results presented in Figs.[5]and[3ta-c), one could 
evaluate the contributions ST?, to T q , within the different 
BZ areas, due to terms including different forms of A q 
in Eq. (IA4|) . The T- and w-dependencies of these 8T q 
terms depend on those of A q , on the integration limits in 
Eq. (|A4I) . and on the number of nonzero integrals there. 
Within the antinodal areas, significant values of ST 9 , are 
obtained through Eq. (|A4I) for uj — » 0, while close to the 
lines of nodes (where T^k) to Eq. (|A3|) is negligible) 
they contribute significantly only when |w| exceeds an 
energy which is somewhat smaller than w A (Q). Thus, 
this energy determines the low-energy scale w L (k) (see 
above) around the lines of nodes. 

As was discussed above, for k B T & \uj\ <c u; L (k'), A q e 
in this k' point is approximately independent of T and 
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U3, within the low-o; A integration sections in Eq. (|A4I) . 
resulting in weak T and u> dependencies of the corre- 
sponding 8T\. For k B T > w L (k') and |w| <C fc B T, since 
Af , self-consistently, scales as 1/T (see above), Jr^j ap- 



proximately scales with T by Eq. (|A4[) . For 3> w L (k') 
and fc B T <C |w|, ^4* e can be approximated in Eq. (|A4p 
by a (^-function around the energy of the corresponding 
arcon or stripon peak, resulting in ST 1 ^ which is linear 
in u>. The high-] uj | extent of this term is limited by the 
high-energy extent of the lagron spectrum in Fig. [5J and 
the high-T extent is limited by the phase stability. 

When k B T or \oj\ exceeds w H , A^ e (due to the 
humpons — see above) behaves similarly to A^ e , and the 
SI 'I term it contributes also has, in this range, a linear 
dependence on T in the low-|a;| limit, and a linear de- 
pendence on u>, in the low-T limit. This results in an 
increase in the slope of the T q vs T (\ui\) curve as the 
increasing k B T (|w|) approaches w H . A similar slope in- 
crease, due to A^ B , is expected when their high-energy 
maxima are approached, as long as they do not exceed 
the high-energy extent of the lagron spectrum. However, 
since (see Fig.^a-b)) there is an apparent discontinuity 
between the QE low- and high-energy BZ areas, where 
the energies at the maxima of A^ ic mostly exceed those 
of the lagron spectrum, the derived contribution of A^ e 
to T q is minute. 



Appendix B: Evaluation of the svivon Scattering 
Rates 

The cosh 2 (££, (w')) and sinh 2 (££, (u/)) factors, appear- 
ing in the expression for T\ in Eq. (f3"T|). are close to 
cosh (2££, (w')) within the LE and ME svivon BZ areas, 
where the major contribution to the expression comes 
from; thus their T and lo dependencies there are close to 
those of cosh (2££,(a/)), derived above through Eq. (131)1) . 
The effect of their oc 1/(T- T ) behavior within the LE 
areas is considered together with the size oc (T — Tq) 2 
of these areas, resulting in an approximate scaling with 
(T — To). The associated decrease in the size of the ME 
areas introduces in these areas a minor oc (T — To) 2 term, 
as in cosh(2££,(a/)). 

The values of n s (k) are determined through Eq. (12U1) . 
and their variation with T takes place mainly in the LE 
BZ areas. However, since n s is T independent, one gets 
through Eq. (|2"2f that the T variation is minor for the av- 
eraged values of n s (k) over the LE areas and the crossover 
areas between them and the ME areas. This applies also 
for the T dependence of the |m s (k')| factors, appearing 
in the expression for T s q in Eq. (13"T1) , due to the proximity 
in the LE areas between the expressions for n s (k) and 
m s (k) in Eq. 

r*(k, uj) represents the scattering of a svivon into a 
particle-hole QE pair, while its spin is condensed; at low 
T and w, it is determined by the stripon and arcon Af e 
spectral functions. In the k B T & \u\ -C regime 
(see Eq. (jA3l) and the discussion in Appendix A), Af c 



is, approximately, constant within the integration range 
in Eq. (13T|) ; in the LE areas, this results in which ap- 
proximately scales with uj/(T — T ), for \u\ < ef nin , its 
magnitude passes through a maximum at higher \u>\ val- 
ues, and crosses over to scaling with w in the \uj\ ^> e^ in 
limit, if it could be reached in this regime; as was dis- 
cussed above, in processes where is integrated over a 
part of the BZ, the contribution of the LE areas to it 
scales with w(T — Tq), for \uj\ < e* lin ; |r®| is smaller, in 
this regime, in the ME areas (and much smaller in the 
HE areas), and approximately scales there with w, plus 
a minor term oc oj(T — Tq) 2 (and it may pass through a 
maximum) . 

When k B T > r^ a & the T-scaled broadening of 
Af c (ui) (see Appendix A) results in a l/T 2 -scaled decrease 
of the above values of T s q with T. For |w| > T q 0a & k B T, 
the contribution of particle- hole QE pairs, based on Af e , 
to is vanishing, resulting in its sharp decrease; a high- 
u) contribution to Y s q due to humpons and high-energy 
QEs is small compared to T s x , for such values of uj. 

The major physical effect of T s q is the introduction, 
right above [4l| Tg, of a width |T S | to the svivon states 
around their energy minima which extends to exceedingly 
low (nonzero) energies \uj\. The corresponding spectral 
functions, obtained through Eq. (|28p. can be approxi- 
mately expressed as: 



A s (k,u) 



r s (k,w)/27r 



[w-e s (k)] 2 + [ir-(k,w)] 2 

r s (k, w ) 



27r[e s (k)] 2 



, for |w| < e s (k). (Bl) 



The oj scaling of T q , at low ui and T, results in an lj- 
independent product A s (k, uj)b T {uj) in the low-w limit. 

Similarly to T q in Eq. (jAl[) . T% could be expressed as 
a sum of two terms: 



n(k,a;)=r;(k,a;)+r$(k,fa;), 



(B2) 



corresponding to the contributions of the Q m lagrons 
(rf), and of the other lagrons (r|), to the q summa- 
tion in Eq. ([3"7|. As in the case of rf in Eq. (fA51) . Tf can 
be approximately expressed as: 



r;(k,w) oc coBh(2fJ(w)) 



(B3) 



x ^[^(k-Q m , w )cosh 2 (^_ Qm ^)) 

m— 1 

-A s (k - Q m , -w) sinh 2 (^_ Qm (-w))]. 

Thus, Tf (k, co)\ has maxima combined of the maxima of 
A s (k — Q TO , |w|) at \uj\ = e s (k - Q m ), for the four Q ?Tl 
points, and decreases to zero at higher energies. Within 
a svivon condensate, e s (k — Q TO ) ~ e s (k) for one of the 
Q m points, and |rj(k,o>)| is largest when both k and 
k — Q m are at the LE areas. In the low uj/k B T limit, Tf 
is negligible (see Eqs. (jBll IB3I0 compared to T q and T| 
(see below). 
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An expression for r|(k, uj) is derived analogously to 
the one for r^k, uj) in Eq. (IA4|) . The q summation 
in Eq. ((3Tj) is, similarly, carried out by using a BZ 
mesh based on circular segments s, where the A s (k — 
q, uj^uj x (q)) cosh 2 (££_ q (w =p <^ A (q)) and A s (k-q, -ui± 
w A (q))sinh 2 (££_ q (-cj ± w A (q)) terms in Eq. d^TJ) are 
expressed as A s s {uj =p w A (q)) cosh 2 (£|(oj =p w A (q)) and 
A s s {— uj ± w A (q)) sinh 2 (£f (-w ± wA (q))> respectively. 

As in the case of T^, the q integral over section s is, 
approximately, proportional to an integral on duj x , where 
the integrand is multiplied by (uj x — w A s s ), between pos- 
itive limits w A | and w A |; again, one has w A s s = uj x ^ = 
within four circular BZ sections around the V-shape la- 
gron energy minima at the Q m points, and |w A s s | are large 
in integration sections corresponding to Q-ESP and high- 
energy SP lagrons. 



Similarly to the evaluation of T\ in Appendix A, the 
b T (uj') functions in Eq. ([57)1 (where uj' is either uj x or 
uj x =F uj) are approximated through their asymptotic be- 
haviors: (a) b T (uj') = k B T/uj, for \uj'\ < k B T; (b) 
b T {uj') S -1, for w' < -fc B T; (c) & r (a/) = 0, for 
a/ 3> & B T- In order to obtain the major features of 
r|(k, uj) in the low- and high-|u;|/£; B T limits, it is for- 
mulated in term of expressions where the asymptotic be- 
havior of b T in (a) is extended to \uj'\ < k B T, and its 
behaviors in (b) and (c) are extended to uj' < —k B T 
and uj' > k B T, respectively. Furthermore, when both 
< uj x < k B T and < uj x T w < k B T, thc differ- 
ences b T (uj x ) — b T (uj x =F uj) are approximated (through 
derivation) as =Fwfc B T '/{uj x =F w /2] 2 - Consequently, one 
gets through Eq. ([37]) : 



rs(k, w ) 



= cosh(2^M)^[r^(k, w )+r^(k, w )+rl(k, w )], where 



duj x {\ - uj^/uj x )[A s s (uj - u/sign(w)) cosh 2 (4 s (w - w A sign(w))) - A^(w A sign(w) - uj) 



xsinh 2 (4 s (o; A sign(w) - uj))} - 



duj x 



[4>-w A sign( w )) 



Ks(Koj) 



x cosh 2 (£ s s (u; - w A sign(w))) - Ag(w A sign(w) - uj) sinh 2 (£ s s (w A sign(w) - uj)) 
du x (l - uJqI/uj x )[A s s {uj + w A sign(a;)) cosh 2 (£f (w + w A sign(a;))) 

-^(-w A sign(w) - uj) sinh 2 (4 s (-w A sign(w) - uj))} 
oc \uj\Kt\- P" duj x {l f - ^ [A^uj - uj x S ign(uj)) cosh 2 (4> - uj x sign(uj))) 

-At (a/signH - uj) sinh 2 (^( W A sign( W ) - uj))} + P" duj x ^ML [A s s (u> + u; A sign(u,)) 



,., (uj x + \uj\/2)* 

x cosh 2 (£ s s (w + w A sign(w))) - A^(-w A sign(w) - uj) sinh 2 (£;(-w A sign(w) \. 

T s us (k,u) (x / ^ A ( W A - Wo As )[A s s (u;-u; A signH)cosh 2 (^( W - W A signH)) 
-^(w A sign(w) - uj) sinh 2 (£f (w A sign(w) - uj))}, and 



(B4) 



UJ 



UJ 



UJ. 



Tl 
°T3 

J T4 

As 
gns 

, Xs 
J 2ns 

Xs 
J 4ns 

Xs 
(ins 



min(|w|,fc B T), uj T2 = max (0, \u\ - k B T), uj a * s = min (uj T2 , oj A | ) , 
max (\oj\,k B T), uj x % = min uj T3 , uj Xs s ) , uj^ is = max (uj x % , w Xs s ) , 

LI i 1. T ,, _™„„/Ti l* T 1 I, ,n , ,As _„;„(',, , ,As 



+ fc B T, uj T5 = max(0,fc B T - |w|), w e * s = min (w T5 ,uO, w /« s = max ( w e4 , w„J) , 
min (w T1 , uj Xs s ) , w, A , s ls = max (uj Xs ns , w A s s ) , uj a ^ s = min (w T1 , w A s s ) , uj a * s = min (w T1 , w A s s ) 
min(|cj|,w A 4), uj^ s = min(|w|,w A s s ), uj x * s = min (w T4 , w^ s ), w A | s = min (w T4 , w A s s ), 



mm (k B T,uj Xs ns ), 



Xs 
fr, 

,Xs\ 



y 4xs 
.As 



min (fc B T, w A s s ), 



J 5ns 
Xs 



mm(k B T, uj^ s ns ), 



J 5xs 



= min (k B T,uj x l), 



mm (w T5 , uj n s s ) , w 6 | s = min ( w T5 , ujJ ) , w 7 * s = min (uj T2 , w„f ) , w 7 | s = min (w T2 , w xs s ) . 
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The dominant contribution to T^, in the low \oj\/k B T 
limit, is of the terms in Eq. (IB4|) . By Fig. [2] their inte- 
gration ranges which contribute significantly to r|(k, w) 
are, for k B T <~w A (Q), within the cone- like areas around 
those Q m points where k — Q TO is within the LE areas; 
the shape of the integration ranges is modified when k B T 
exceeds ~cj a (Q). 

As was mentioned above, at a point k' within the LE 
areas, A s (k',uj') is not small in the range w' ~ k B T, due 
the effect of F* at T > T«; by Eq. O, its T dependence 
could then be approximated (within the range where it is 
significant) in terms of a scaling factor a s between u' and 
(T — To) which yields, taking into account normalization: 



(T' - T )A s (k',uj' = ±a s (T' - T )) I 
s (T" -T )A s (k',u>' = ±a s (T" -T )) 
for k' G the LE areas. 



T = T' 
a T = T", 

(B5) 



Since the cosh 2 and sinh 2 factors in the expression for 
r* s , in Eq. (|B4|) . approximately scale in the LE areas with 
1/(T-T ), they introduce a 1/(T-T ) factor to the in- 
tegrals within these areas. Furthermore, since the size of 
the LE areas scales with (T-T ) 2 , and (T-T )A s (k', w A ) 
has a scaling dependence on cj a / (T — To) (see Eq. (|B5[) ). 
one could, approximately, replace the variable lo x in these 
integrals by ui x /(T — Tq), in the low \oj\/k B T limit at least 
for k B T <-w A (Q). 

This results in r|(k, uS) which approximately scales, in 
the low \u>\/k B T limit, with w/[(T - T ) 2 (l - T /T)], if k 
is within the LE areas, and with cj/[(T - T )(l - T /T)}, 
if it is within the part of ME areas where there exists a 
Q m point for which k — Q m is within the LE areas (which 
applies when k B T <^w A (Q)). Thus, the u dependence 
of r|(k, oj) in the low \u>\/k B T limit is similar to that of 
r^(k,w) (see above) in the k B T <C Tq & regime, but it 
decreases faster with T there; however, this behavior of 
r| persists in the k B T » regime, where decreases 
faster with T, and it is somewhat modified when k B T 
exceeds ~w A (Q). 

Similarly to T s q above, an integration on the low- 
(|w|/fc B r) r|, over a part of the BZ, results in a con- 
tribution of the LE areas which, approximately, scales 
with w/(l - T /T): for k B T >~w A (Q) this contribution 
persists as the dominant part of the k-integrated low- 
(\uj\/k B T) r s , and it approximately scales with u. 

In the high-|w|/A: B T limit, T| is determined by the r* s 
terms in Eq. (|B4[) . The A s functions in the integrals there 
can then be approximated by (5-functions around the cor- 
responding ±e s (k'), resulting in factors \ui\ — e s (k') — Wq s s 
(thus depending linearly on ui), for |w| > e s (k'); however, 
there are also multiplicative factors of cosh 2 (££, (e s (k')) 
or sinh 2 (££, (e s (k')) and of cosh (2££(w)) which has an 
opposite w-dependence for \u>\ <^J. 



Appendix C: Evaluation of the Electron Spectral 
Functions 

By Eq. (|4"2"|) , Ag(k, cj) is a convolution of QE and svivon 
spectral functions, weighted by the b T (uj') + f T (ui' =F w ) 
factors. In order to study the effect of these factors, 
let us approximate them through their low- and high- 
\uj'\/k B T asymptotic behaviors; thus they are expressed 
as: (a) b T (a/) + f T (u>' T w) = k B T/u}', for < k B T; (b) 
b T {uj') + / T (a/ T w) = 0, for |w'| & |a/ T w| > k B T and 
sign(w') = sign(o;'=Fw); (c) 6 t (w') + /tKt u ) = k//V|, 
for \ui'\ & |u/ =F uj\ > k B T and sign (a/) ^ sign(w' =F u); 
(d) b T {oj') + / T (u/ T w) = o//2|o/|, for > k B T and 
|u/ T w| < k B T. Eq. (|42]) is then approximated as: 

A d (k,cj) = 4 b (k,c)+<(k, W ), where 

^ k' W 

X [A 9 (k - k', CJ - U)') cosh 2 (^,(w')) 
+A 9 (k - k', u + J) sinh 2 (& (w'))] , 

<(k, w ) S -£<| / doA4«(k>') 



x[l-^(* B T-H+u')] 

x [6»(w)A 9 (k -k',w - u/)cosh 2 (££,(w')) 

+6»(-o;) J 4«(k - k', u + w') sinh 2 (££, (a/))] 



(CI) 



du/A s (k',c/) 

'-M-fc B T 

xfi-^r-H-w')] 

x [0(-w)A»(k - k', w - u/) cosh 2 (^,(w')) 
+6(uj)A c '(k -X,u + u/) sinh 2 (££,(u/))] 

^o c (k, ui) in Eq. (IC1|) represents a convolution of QE 
and modified svivon spectral functions, where the low- 
energy (\uj'\ < k B T) tails of A s have been truncated (see 
Eq. (|B1I) and Fig. [Bp})), while A^ h (k,uj) is generated by 
these A s tails and QE spectral functions. Since the QE 
spectrum (see Fig.[3]Ja-c)) includes almost dispersionless 
low-energy peaks extending over ranges of the BZ, the 
effect of the A q (k — k',ui =p u/) terms in the expression 
for A^ h is to introduce their low-energy peaks which are 
modified (due to the k' summation and oj' integration) 
to increase the peak widths by ~2fc B T; this results in 
an effective low-energy (LE) electron band of a linewidth 
which depends linearly on T, for T > 

The contribution of svivons to the k' summation (in- 
cluding averaging over the svivon condensates) in the ex- 
pression for A^ h in Eq. (|C1I) is determined by that of the 
LE svivons, plus some contribution of ME svivons. For 
LE svivons, the cosh 2 and sinh 2 factors in this expression 
approximately scale with 1/(T — Tq) (see Eq. ([391) ). and 
the T dependence of the A s (k' , ui') terms there is approx- 
imated through the scaling relation specified in Eq. (|B5I) , 
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introducing a factor of 1/(T — To), and replacing the in- 
tegration variable uj' by x — us' /k B {T — T ) (with inte- 
gration limits ±T/(T - T )). For ME svivons the cosh 2 , 
sinh 2 and A 5 factors have a weak T dependence. 

As was discussed above, the size of the BZ areas of 
LE svivons scales with (T — To) 2 (on the expense of ar- 
eas of the ME svivons). The expression for ^Q b (k, uj) in 
Eq. (|C1[) introduces, for k' points corresponding to LE 
svivons, a BZ section (of the same size) of k — k' points of 
QEs, out of which only its subsection of low-energy QEs 
contributes to the LE A°\ h effective band. The size of 
this subsection approximately scales with (T — To)" 9 ^, 
where the T-dependent exponent a q (k) > is generally 
< 2. Consequently, the k' summation over LE svivons 
in Eq. jCT|) introduces a factor of (T - To) a " {ii) to their 
contribution to the LE Ap b (k, ui). 

At low T > T q , these low-QE-energy k — k' points are 
shifted from the electron k point (choosing the QE BZ 
as in Fig. |^a-c)) by ~ (<5q m /2) and ~ (Q m - <5q m /2), 
for the four values of m (see Eq. (|2T)1) ). As can be viewed 
in Fig. [3jc) , a q (k) is close to 2 for antinodal k points, 
within this temperature regime, and one can then also as- 
sume that k B T rp^k), where the T dependence of the 
stripon A q is weak (see Eq. (IA3|) and discussion in Ap- 



pendix A). Thus, one gets by Eq. (|C1|) and the above dis- 
cussion that the antinodal LE Ap b (k, w) approximately 
scales with T around its maximum; since the width of 
the effective electron band, it represents, is characterized 
by a constant plus a linear T term, this implies that for 
antinodal electrons at low T > T q , the spectral weight 
within the LE effective band anomalously increases with 
T. There is some increase with T, in this regime, also in 
the spectral weight within the nodal LE effective band, 
mainly due to the contribution of ME svivons. 

On the other hand, for antinodal and nodal k points 
in the k B T ^> r^(k) regime (in the nodal case this cor- 
responds to T > T q ), a q (k) is close to 1 (see Fig. ^c)), 
and A q in Eq. (|C1[) scales with 1/T (see discussion in Ap- 
pendix A). This results in LE ^Q b (k, ui) which approxi- 
mately scales with 1/T around its maximum; thus since 
the linewidth of the effective band it represents scales 
with T, one gets that the spectral weight within this band 
is approximately T independent for k B T ^> T^k). Note, 
however, that when k B T approaches w H , the humpons 
merge with the stripons in the antinodal BZ areas, result- 
ing in a further increase in the spectral weight within the 
effective LE electron band there, saturating at k B T > cj h . 
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